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Abstract. The availability of the Hipparcos Catalogue triggered many kinematic and dynamical studies of the solar neigh- 
bourhood. Nevertheless, those studies generally lacked the third component of the space velocities, i.e., the radial velocities. 
This work presents the kinematic analysis of 5952 K and 739 M giants in the solar neighbourhood which includes for the first 
time radial velocity data from a large survey performed with the CORAVEL spectrovelocimeter. It also uses proper motions 
from the Tycho-2 catalogue, which are expected to be more accurate than the Hipparcos ones. An important by-product of this 
study is the observed fraction of only 5.7% of spectroscopic binaries among M giants as compared to 13.7% for K giants. After 
excluding the binaries for which no center-of-mass velocity could be estimated, 5311 K and 719 M giants remain in the final 
sample. 

The l/V-plane constructed from these data for the stars with precise parallaxes ((T„/n < 20%) reveals a rich small-scale 
structure, with several clumps corresponding to the Hercules stream, the Sirius moving group, and the Hyades and Pleiades 
superclusters. A maximum-likelihood method, based on a bayesian approach, has been applied to the data, in order to make full 
use of all the available stars (not only those with precise parallaxes) and to derive the kinematic properties of these subgroups. 
Isochrones in the Hertzsprung-Russell diagram reveal a very wide range of ages for stars belonging to these groups. These 
groups are most probably related to the dynamical perturbation by transient spiral waves (as recently modelled by De Simone 
et al. 2004) rather than to cluster remnants. A possible explanation for the presence of young group/clusters in the same area of 
the t/K-plane is that they have been put there by the spiral wave associated with their formation, while the kinematics of the 
older stars of our sample has also been disturbed by the same wave. The emerging picture is thus one of dynamical streams 
pervading the solar neighbourhood and travelling in the Galaxy with similar space velocities. The term dynamical stream is 
more appropriate than the traditional term supercluster since it involves stars of different ages, not bom at the same place nor at 
the same time. The position of those streams in the L'V-plane is responsible for the vertex deviation of 16.2° ±5.6° for the whole 
sample. Our study suggests that the vertex deviation for younger populations could in fact have the same dynamical origin. 
The underlying velocity ellipsoid, extracted by the maximum-likelihood method after removal of the streams, is not centered 
on the value commonly accepted for the radial antisolar motion: it is centered on {U) = -2.78 ± 1.07 km s"'. However, the 
full data set (including the various streams) does yield the usual value for the radial solar motion, when properly accounting 
for the biases inherent to this kind of analysis (namely, (U) = -10.25 ± 0.15 km s"'). This discrepancy clearly raises the 
essential question of how to derive the solar motion in the presence of dynamical perturbations altering the kinematics of the 
solar neighbourhood: does there exist in the solar neighbourhood a subset of stars having no net radial motion which can be 
used as a reference against which to measure the solar motion? 

Key words. Galaxy: kinematics and dynamics - disk - solar neighbourhood - evolution - structure - stars: kinematics - 
late-type 
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1. Introduction 

The kinematics of stars in the solar neighbourhood, studied 
from proper motions and radial velocities, has long been ac- 
knowledged to be an essential element for the understanding 
of the structure and evolution of the Galaxy as a whole (e.g. 
Schwarzschild 1907, Lindblad 1925, Oort 1927). The story of 
the efforts to obtain the necessary data is standard textbook 
fare, and nicely illustrates how theoretical progress and data 
acquisition have to go hand in hand if one wants to gain in- 
sight in the true nature of our Galaxy. Clearly, another chap- 
ter in this story has begun with the Hipparcos satelUte mission 
(ESA 1997), that provided accurate parallaxes and proper mo- 
tions for a large number of stars (about 118 000): with accu- 
rate parallaxes available, it was no longer necessary to resort 
to photometric distances. The Hipparcos data enabled several 
kinematic studies of the solar neighbourhood, but those stud- 
ies generally lacked radial velocity data. A first preview of how 
Hipparcos data could improve our knowledge of stellar motions 
in the Galaxy was given by Kovalevsky (1998); then, Dehnen & 
Binney (1998) used Hipparcos proper motions to derive some 
fundamental kinematic parameters of the solar neighbourhood. 
They did not use radial velocities from the literature because, at 
that time, radial velocities had been measured preferentially for 
high proper motions stars and their use would have introduced 
a kinematic bias. Hipparcos data were also used by Chereul 
et al. (1998) to derive the small scale structure of the velocity 
distribution of early-type stars in the solar neighbourhood. 

The situation has dramatically improved thanks to the ef- 
forts of a large European consortium to obtain radial velocities 
of Hipparcos stars later than about F5 (Udry et al. 1997). The 
sample includes Hipparcos 'survey' stars (flag S in field H68 
of the Hipparcos Catalogue), and stars from other specific pro- 
grammes. This unique database, comprising about 45 000 stars 
measured with the CORAVEL spectrovelocimeter (Baranne et 
al. 1979) at a typical accuracy of 0.3 km s"\ combines a 
high precision and the absence of kinematic bias. It thus rep- 
resents an unprecedented data set to test the results obtained 
by previous kinematic studies based solely on Hipparcos data. 
Another feature of the present study is the use of Tycho-2 
proper motions, which combine Hipparcos positions with po- 
sitions from much older catalogues (H0g et al. 2000). They 
represent a subsantial improvement over the Hipparcos proper 
motions themselves, which are based on very accurate posi- 
tions, but extending over a limited 3-year time span. Results 
from the Geneva-Copenhagen survey for about 14 000 F and G 
dwarfs present in the CORAVEL database have been recently 
pubUshed by Nordstrom et al. (2004). The present paper, on 
the other hand, concentrates on the K and M giants from the 
CORAVEL database. 

According to the classical theory of moving groups, early- 
type stars are not good probes to determine global character- 
istics of the Galaxy because these stars are young and still 
carry the kinematic signature of their place of birth. As a con- 
sequence, their distribution in velocity space is clumpy (e.g. 
de Bruijne et al. 1997, Figueras et al. 1997, de Zeeuw et al. 
1999). These inhomogeneities can be spatially confined groups 
of young stars (OB associations in the Gould's belt, young 



clusters) but can also be spatially extended groups. There is, 
however, some confusion in the literature about the related ter- 
minology. Eggen (1994) defines a 'supercluster' as a group 
of stars gravitationally unbound that share the same kinemat- 
ics and may occupy extended regions in the Galaxy, and a 
'moving group' as the part of the supercluster that enters the 
solar neighbourhood and can be observed all over the sky. 
Unfortunately, the same term 'moving group' is sometimes 
also applied to OB associations (e.g. de Zeeuw et al. 1999). 
It has long been known that, in the solar vicinity, there are 
several superclusters and moving groups that share the same 
space motions as well-known open clusters (Eggen 1958). The 
best documented groups (see Montes et al. 2001 and references 
therein) are the Hyades supercluster associated with the Hyades 
cluster (600 Myr) and the Ursa Major group (also known as 
the Sirius supercluster) associated with the UMa cluster of 
stars (300 Myr). Another kinematic group called the Local 
Association or Pleiades moving group is a reasonably coher- 
ent kinematic group with embedded young clusters and asso- 
ciations like the Pleiades, a Persei, NGC 2516, IC 2602 and 
Scorpius-Centaurus, with ages ranging from about 20 to 150 
Myr. Two other young moving groups are the IC 2391 super- 
cluster (35-55 Myr) and the Castor moving group (200 Myr). 
The kinematic properties of all these moving groups and su- 
perclusters are listed by Montes et al. (2001) and Chereul et 
al. (1999). However, the recent observation (e.g. Chereul et al. 
1998, 1999) that some of those superclusters involve early-type 
stars spanning a wide range of ages contradicts the classical hy- 
pothesis that supercluster stars share a common origin: Chereul 
et al. (1998) propose that the supercluster- like velocity struc- 
ture is just a chance juxtaposition of several cluster remnants. 

Nevertheless, significant dumpiness in velocity space has 
also been recently reported for late-type stars (Dehnen 1998), 
thus raising the questions of the age of those late-type stars and 
of the exact origin of this dumpiness. For example, Montes et 
al. (2001) suppose that, if a late-type star belongs to a moving 
group or supercluster, it may be considered as a sign that the 
star is young, an assumption that will be largely challenged by 
the results of this paper. Here, we attack the problem from an- 
other point of view by restricting our study of stellar kinematics 
to late-type giants (of spectral types K and M), without any a 
priori hypothesis on their motion or age. 

The present paper is organized as follows. Sect. 2 describes 
the various data sets used in the present study. The kinematic 
analysis of the data and its implications for the dynamics of our 
Galaxy are presented in Sect. 3. For summary and conclusions 
we refer to Sect. 4. 

2. The data 

2.1. Selection criteria 

Our stellar sample is the intersection of several data sets: (i) in 
a first step, all stars with spectral types K and M appearing in 
field H76 of the Hipparcos Catalogue (ESA 1997) have been 
selected, and the Hipparcos parallaxes have been used in the 
present study; (ii) Proper motions were taken from the Tycho-2 
catalogue (H0g et al. 2000). These proper motions are more ac- 
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curate than the Hipparcos ones. Fig.ncompares the proper mo- 
tions from Tycho-2 and Hipparcos, and reveals that the Tycho- 
2 proper motions are sufficiently different from the Hipparcos 
ones (especially for binaries) to warrant the use of the former 
in the present study; (iii) Radial-velocity data for stars belong- 
ing to this first list (stars with spectral type K and M) have 
then been retrieved from the CORAVEL database. Only stars 
from the northern hemisphere (6 > 0°), observed with the 1-m 
Swiss telescope at the Haute-Provence Observatory, have been 
considered. As described by Udry et al. (1997), all stars from 
the Hipparcos survey (including all stars brighter than V - 
13 + 1.1| sin^l for spectral types later than G5, where b is the 
galactic latitude; those stars are flagged 'S' in field H68 of the 
Hipparcos Catalogue) are present in the CORAVEL database. 
Besides the survey, stars from several other Hipparcos pro- 
grammes, dealing with galactic studies and requiring radial ve- 
locities, were also monitored with CORAVEL. The CORAVEL 
database finally contains stars monitored for other purposes 
(e.g., binarity or rotation), but they represent a small fraction 
of the Hipparcos survey stars. 

The sample resulting from the intersection of these 3 data 
sets was further cleaned as follows. First, a crude Hertzsprung- 
Russell diagram for all the M-type stars of the Hipparcos 
Catalogue has been constructed from distances estimated from 
a simple inversion of the parallax (Fig. IS, and reveals a clear 
separation between dwarfs and giants. Clearly, all M stars with 
Mhp < 4 must be giants. For K-type stars with a relative par- 
allax error less than 20%, the Hertzsprung-Russell diagram 
(Fig. O shows that the giant branch connects to the main se- 
quence around {V - I - 0.75, Mnp - 4). Therefore, in order 
to select only K and M giants, only stars with Mnp < 4 were 
kept in the sample. For K stars, the stars with V - I < 0.75 and 
Mhp > 2 were eliminated to avoid contamination by K dwarfs. 
Diagnostics based on the radial-velocity variability and on the 
CORAVEL cross-correlation profiles, combined with literature 
searches, allowed us to further screen out stars with peculiar 
spectra, such as T Tau stars, Mira variables and S stars. T Tau 
stars have been eliminated because they belong to a specific 
young population; S stars because they are a mixture of extrin- 
sic and intrinsic stars, which have different population char- 
acteristics (Van Eck & Jorissen 2000); Mira variables because 
their center-of-mass radial velocity is difficult to derive from 
the optical spectrum, where confusion is introduced by the pul- 
sation (Alvarez et al. 2001). 

The primary sample includes 5952 K giants and 739 M 
giants (6691 stars). 86% of those stars are 'survey' stars: our 
sample is thus complete for the K and M giants brighter than 
V - 7.3-H 1.1| sin/?!. To fix the ideas, for a typical giant star with 
My - 0, this magnitude threshold translates into distances of 
290 pc in the galactic plane and 480 pc in the direction of the 
galactic pole (these values are, however, very sensitive upon 
the adopted absolute magnitude, and become 45 pc for a sub- 
giant star with My - +4 in the galactic plane, and 2900 pc for 
a supergiant star with My - -5, corresponding to the range 
of luminosities present in our sample; see Sect. I3.3.6l l. The fi- 
nal step in the preparation and cleaning of the sample involves 
the identification of the spectroscopic binaries, as described in 
Sect.lO 
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Fig. 1. Difference between the Hipparcos and Tycho-2 proper- 
motion moduli, normalized by the root mean square of the stan- 
dard errors on the Hipparcos and Tycho-2 proper motions, de- 
noted e^. The solid line refers to the 859 spectroscopic bina- 
ries (SB) present in our sample, whereas the dashed line cor- 
responds to the 5832 non-SB stars. Note that the sample of bi- 
naries contains more cases where the Hipparcos and Tycho-2 
proper motions differ significantly, in agreement with the argu- 
ment of Kaplan & Makarov (2003). 
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Fig. 2. Crude Hertzsprung-Russell diagram for the Hipparcos 
M stars with positive parallaxes. 

2.2. Binaries 

The identification of the binaries, especially those with large 
velocity amplitudes, is an important step in the selection pro- 
cess, because kinematic studies should make use of the center- 
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Fig. 3. Crude Hertzsprung-Russell diagram for the Hipparcos 
K stars with a relative error on the parallax less than 20%. 



of-mass velocity. In order to identify those, the observing strat- 
egy was to obtain at least two radial- velocity measurements per 
star, spanning 2 to 3 yr. Monte-Carlo simulations reveal that, 
with such a strategy, binaries are detected with an efficiency 
better than 50 percent (Udry et al. 1997). 

Since late-type giants exhibit intrinsic radial velocity jitter 
(Van Eck & Jorissen 2000), the identification of binaries re- 
quires a specific strategy, which makes use of the {S b, croivro)) 
diagram (Fig.|5}, where vvo is the measured radial velocity. The 
parameter Sb is a measure of the intrinsic width of the cross- 
correlation profile, i.e., corrected from the instrumental width 
(see Fig.l^for more precise definitions). The CORAVEL instru- 
mental profile is obtained from the observation of the cross- 
correlation dip of minor planets reflecting the sun light, after 
correction for the solar rotational velocity and photospheric tur- 
bulence (see Benz & Mayor 1981 for more details). The Sb 
parameter is directly related to the average spectral line width, 
which is in turn a function of spectral type and luminosity, the 
later-type and more luminous stars having larger 5 values. On 
the other hand, the measurement error e has been quadratically 
subtracted from the radial velocity standard deviation cr(v,o) to 
yield the eff'ective standard deviation cro(vro). 

The identification of binaries among K or M giants fol- 
lows diff'erent steps. The standard variability test, compar- 
ing the standard deviation cr(vro) of the measurements to their 
uncertainty e, cannot be applied to M giants, because intrin- 
sic radial-velocity variations ('jitter') associated with envelope 
pulsations would flag them as velocity variables in almost all 
cases. Therefore, in a first step, M giants having croivro) > 
1 km s 'have been monitored with the ELODIE spectro- 
graph (Baranne et al. 1996) at the Haute-Provence Observatory 
(France) since August 2000 (Jorissen et al. 2004). These sup- 
plementary data points made it possible to distinguish orbital 
variations from radial- velocity jitter by a simple visual inspec- 



tion of the data. It turns out that, in the {S b, cro(vvo)) diagram, 
all the confirmed binaries (filled symbols in the lower panel 
of Fig. |5} are located in the upper left corner, and are clearly 
separated from the bulk of the sample. Stars located below 
the dashed line in Fig. [S] may be supposed to be single (al- 
though some very long -period binaries, with P t 10 y, may 
still hide among those). Their radial velocity dispersion suf- 
fers from a jitter which clearly increases with increasing spec- 
tral line- width, as represented by the parameter Sb. For very 
large values of Sb (in excess of about 9 km s"'), the dia- 
gram is populated almost exclusively by supergiants (star sym- 
bols in the lower panel of Fig. |5}. Many semi-regular and ir- 
regular variables are located in the intermediate region, with 
5 ^ Sb ^ 9 km s"'. Strangely enough, spectroscopic binaries 
seem to be lacking in this region. A more detailed discussion 
of the properties of the binaries found among M giants is de- 
ferred to a forthcoming paper (see Jorissen et al., 2004 for a 
prehminary account). 

For K giants, no such structure is apparent in the 
(S b,o-()(vro)) diagram. It has been checked that the distribu- 
tion of stars along a (N - l)[o"(v',o)/e]- axis (where is the 
number of measurements for a given star) follows a distri- 
bution, as expected (e.g., Jorissen & Mayor 1988). This holds 
true irrespective of the S b value, thus confirming the absence 
of structure in the (S b, cro(vro)) diagram. The binaries among K 
giants may thus be identified by a straight test. A 1 % thresh- 
old for the first kind risk (of rejecting the null hypothesis that 
the star is not a binary while actually true, i.e., of considering 
the star as binary while actually single) has been chosen in the 
present study. 

Among M giants, 42 binaries are found (corresponding to 
an observed frequency of spectroscopic binaries of 42/739, or 
5.7%). Among the 5952 K giants, 817 spectroscopic binaries 
are found, corresponding to a frequency of 13.7%. The large 
diff'erence between these two frequencies hints at some under- 
lying physical cause, to be discussed in a forthcoming paper 
(see Jorissen et al. 2004 for a preliminary discussion). Most of 
the spectroscopic binaries (SBs) detected among our samples 
of K and M giants are first detections. This large list of new 
SBs constitutes an important by-product of the present work. 
The new binaries are identified in Table A.l (flags 0, 1, 5, 6 
and 9 in column 24 of Table A.l) and will be the topic of a 
separate paper 

Among this total sample of 859 spectroscopic binaries, 
the center-of-mass velocity could be computed (whenever the 
available measurements were numerous enough to derive an 
orbit), estimated (in the case of low-amplitude orbits) or taken 
from the literature for 216 systems only, thus leaving 643 bi- 
naries which had to be discarded from the kinematic study, be- 
cause no reliable center-of-mass velocity could be estimated. 

After excluding those large-amplitude binaries as well as 
the dubious cases, 5311 K giants and 719 M giants remain in 
the final sample. Fig.|6lshows the distribution of our final sam- 
ple on the sky. Among this final sample of 6030 stars, 5397 
belong to the Hipparcos 'survey'. 
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Fig. 5. The (Sb,o-o(vro)) diagram (see text) for K (upper panel, and zoom inside) and M giants (lower panel). Star symbols 
denote supergiant stars, filled symbols denote spectroscopic binaries with (squares) or without (circles) center-of-mass velocity 
available. 



3. Kinematic analysis 

In this section, we present the kinematic analysis of the data. 
Although the velocity distribution is a function of position in 
the Galaxy, it will be assumed here that our local stellar sam- 
ple may be used to derive the properties of the velocity distri- 
bution at the location of the Sun if we correct the data from 
the effects of differential galactic rotation. First, we analyze 
the kinematics of the sample restricted to the 2774 stars with 
parallaxes accurate to better than 20%. In a second approach, 
designed to make full use of the 6030 available stars but with- 
out being affected by the biases appearing when deaUng with 
low-precision parallaxes, the kinematic parameters are evalu- 
ated with a Monte-Carlo method. Although easy to implement, 
the Monte Carlo method faces some Umitations (the parallaxes 
were drawn from a Gaussian distribution centered on the ob- 
served parallax, not the true one as it should). Therefore, in a 



third approach, we make use of a bayesian method (the LM 
method; Luri et al. 1996), which allows us to derive simul- 
taneously maximum likelihood estimators of luminosity and 
kinematic parameters, and which can identify possible groups 
present in the sample by performing a cluster analysis. 

In all cases, we correct for differential Galactic rotation in 
order to deal with the true velocity of a star relative to the Local 
Standard of Rest (the LSR is a reference frame in circular orbit 
in the Galactic plane at the galactocentric distance of the Sun). 
Given the observed values of the radial velocity Vro and the 
proper motions jU;o and fibo, the corrected values are (Trumpler 
& Weaver 1953): 

Vr = Vro -Adcos^bsitiZl (1) 
K Hi = Kfiio - A cos2l - B (2) 
K/^b = KHbo + 1/2 A sin2fo sin2/, (3) 
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Fig. 4. The concept of line-width parameter 5 = {s^ - s^y^^ is 
illustrated here by comparing the CORAVEL cross-correlation 
(smoothed) profiles (black dots) of 3 M giants or supergiants 
with the gaussian instrumental profile (of sigma so -7 km s"' ; 
solid line). Each profile coiTesponds to a single radial-velocity 
measurement for the given stars. The s value listed in the above 
panels refers to the sigma parameter of a gaussian fitted to the 
observed coiTelation profile. 




Fig. 6. Distribution of the final sample on the sky, in galactic 
coordinates, the galactic center being at the center of the map. 
The selection criterion 5 > 0° is clearly apparent on this map. 



where k is the factor to convert proper motions into space ve- 
locities, and d is the distance to the Sun. We used the values 
of the Oort constants A and B derived by Feast & Whitelock 
(1997) from Hipparcos Cepheids, i.e. A = 14.82 kms^'kpc ' 
and B = -12.37 kms 'kpc"'. The observed velocity is, after 
that correction, a combination of the peculiar velocity of the 
star and that of the Sun relative to the LSR. 

3.1. Analysis of the sample restricted to stars with the 
most precise parallaxes 

The components of the velocity of a star with respect to the 
Sun in the cartesian coordinate system associated with galactic 
coordinates are the velocity towards the galactic center U, the 
velocity in the direction of Galactic rotation V, and the verti- 



cal velocity W. The basic technique to calculate U, V and W is 
to invert the parallax to estimate the distance and then use the 
proper motions and radial velocities as given by Eqs. (1), (2) 
and (3). This simple procedure faces two major difficulties: (i) 
the inverse parallax is a biased estimator of the distance (espe- 
cially when the relative error on the parallax is high), and (ii) 
the individual errors on the velocities cannot be derived from a 
simple first-order linear propagation of the individual errors on 
the parallaxes. It is thus extremely hard to estimate an error on 
the average velocities. If we limit the sample to the stars with 
a relative parallax eiTor smaller than 10%, we are left with 786 
stars, which is too small a sample to analyze the general be- 
haviour of the K and M giants in the solar neighbourhood. We 
choose instead to restrict the sample to the stars with a relative 
parallax error smaller than 20%, in such a way that 2774 stars 
(2524 K and 250 M) remain. In that case the classical first order 
approximation for the calculation of the errors may still be ap- 
plied (and the bias is very small, see Brown et al. 1997), which 
yields from the errors on the proper motions, radial velocities 
and parallaxes 

(eu) =4.03kms-i 

(ev) = 3.22kms-i (4) 
(ew) = 2.54 km s"'. 

The contribution of the measurement eiTors to the uncer- 
tainty on the sample mean velocity is N^^^^ times the values 
given by ©, where A^(= 211 A) is the sample size. This contri- 
bution is in fact negligible with respect to the "Poisson noise" 
(obtained from the intrinsic velocity dispersion of the sample, 
see Eq. (|8}). For the stars with a relative parallax error smaller 
than 20%, we obtain 

{U) = -10.24 + 0.66 kms-i 

(V) = -20.51 +0.43 km s-i (5) 
{W) = -7.77 + 0.34 km s-i. 

If there is no net radial and vertical motion at the solar po- 
sition in the Galaxy, we have hence estimated Uq - -{U) and 
Wo = -{W) (see however Sect. 3.3.6). 

Knowing that (V) is affected by the asymmetric drift, which 
implies that the larger a stellar sample's velocity dispersion is, 
the more it lags behind the circular Galactic rotation, it is inter- 
esting to compare the mean value of V for the K and M giants 
separately. For the 250 M giants, we obtain 

{V) = -23.42 + 1.48 km s"' (6) 

while for the 2524 K giants we obtain 

{V) = -20.22 + 0.44 km s"'. (7) 

This difference (at the 2cr level) between the two subsamples 
can be understood in terms of the age-velocity dispersion re- 
lation. Indeed, the M giants must be a Uttle older than the K 
giants on average because only the low-mass stars can reach 
the spectral type M on the Red Giant Branch and because the 
lifetime on the main sequence is longer for lower mass stars. 
This implies that the subsample of M giants has a larger ve- 
locity dispersion and rotates more slowly about the Galactic 
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Center than the subsample of K giants, in agreement with the 
asymmetric drift relation. 

The velocity dispersion tensor is defined as 
((v - (v)) <8i (v - <v))). The diagonal components are the 
square of the velocity dispersions while the mixed components 
correspond to the covariances. For the stars with a relative 
parallax error smaller than 20%, we obtain the classical 
ordering for the diagonal components, i.e. (Ju > a-y > cw'- 



o-u = 34.46 + 0.46 km S-' 
o-v = 22.54 ±0.30 km s-i 
o-w = 17.96 + 0.24 km s-i. 



(8) 



The asymmetric drift relation predicts a linear dependence of 
(v) - {V} + Vq with the radial velocity dispersion cr^. If we 
adopt the peculiar velocity of the sun from Dehnen & Binney 
(1998), Ve = 5.25 kms ', we find for the full sample (K and M 
together) 



(v) = -15.26kms" 



(9) 



On the other hand, if we adopt the parameter ^ = 80 + 5 km s ' 
from the asymmetric drift equation of Dehnen & Binney 
(1998), we find 

(v) = -o-^/;t = -14.9 + 1.3 kms-i. (10) 

This independent estimate of (v) is in accordance with (|9} and 
our values of (V) and cr^ are thus in good agreement with this 
value of k. 

As we already noticed, the M giants are a little bit older 
than the K giants on average and we thus expect the velocity 
dispersions of the subsample of M giants to be higher We find 



cTf/ (M giants) - 35.95 ± 1.6kms 



and 



(11) 
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Fig. 7. Density of stars with precise parallaxes {cr„/n < 20%) 
in the f/V-plane. The colours indicate the number of stars in 
each bin. The contours indicate the bins with 3, 4, 7, 12, 17, 20, 
25, 30, 35 and 40 stars respectively. The concentration of stars 
around U ^ -35 km s"', V ^ -45 km s"' contributes largely 
to the vertex deviation, while other peaks already identified by 
Dehnen (1998) at t/ = -40 km s"', V = -25 km s ' (Hyades 
supercluster), atU - 10 km s"', V = -5 km s"' (Sirius moving 
group), and at f/ = -15 km s"', V - -25 km s"' (Pleiades 
supercluster) are also present. 



o-u(K giants) = 34.32 ± 0.48 km s" 



As expected, the difference is not very large, and the radial ve- 
locity dispersions for the two subsamples satisfy the asymmet- 
ric drift equation ( I10> . 

The mixed components of the velocity dispersion ten- 
sor involving vertical motions vanish within their errors. 
Nevertheless, the mixed component in the plane, that we de- 
note cr^y, is non-zero: this is not allowed in an axisymmetric 
Galaxy and is further discussed in Sect. 3.3. We obtain 



aly = 134.26+ 13.28 km^s"^ 



where the eiTor coiTesponds to the 15 and 85 percentiles of the 
correlation coeflicient, assuming that the sample of U, V veloc- 
ities is drawn from a two-dimensional gaussian distribution. 

In order to parametrize the deviation from dynamical ax- 
isymmetry, a useful quantity is the vertex deviation /,„ i.e., the 
angle one has to rotate the (U, V) coordinate system in order 
to diagonalize the velocity dispersion tensor The value of the 
vertex deviation for our sample is 



•4)) 



10.85° + 1.62° 



(12) This vertex deviation for giant stars is not in accordance with 
a perfectly axisymmetric Galaxy and could be caused by non- 
axisymmetric perturbations in the solar neighbourhood, or by a 
deviation from equilibrium (i.e. moving groups due to inhomo- 
geneous star formation, following the classical theory of mov- 
ing groups). A hint to the true nature of this vertex deviation 
could be the local anomaly in the f/V-plane, the so-called u- 
anomaly (e.g. Raboud et al. 1998). If we calculate the mean 
velocity (U) of the stars with V < -35kms"', we find that 
it is largely negative ({U} = -22kms"'). It denotes a global 
outward radial motion of the stars that lag behind the galactic 
rotation. We see on Fig.0that this anomaly is due to a clump 
located at U - -35 km s"', V ^ -45 km s"', the already 
known 'Hercules' stream (Fux 2001). On the other hand, Fig.0 
reveals in fact a rich small-scale structure in the f/V-plane, 
with several clumps which can be associated with known kine- 
matic features: we clearly see small peaks at f/ = -40 km s"', 
V = -25 km s"' (corresponding to the Hyades supercluster), 
at U - 10 km s"', V = -5 km s"' (Sirius moving group), and 
at U - -15 km s"', V = -25 km s"' (Pleiades supercluster). 
More precise values for these peaks will be given in Sect. 3.3 
(14) where their origin will also be discussed. 



(13) 
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3.2. Monte Carlo simulation 

To reduce the bias introduced by the non-linearity of the 
parallax-distance transformation, the results discussed in the 
previous Section were obtained by restricting the sample to 
stars with a relative parallax error smaller than 20%. In this 
Section, we use all the available stars, without any parallax 
truncation, but perform a Monte-Carlo simulation to properly 
evaluate the errors and biases on the kinematic parameters. 

We constructed synthetic samples by drawing the stellar 
parallax from a gaussian distribution centered on the observed 
parallax value for the corresponding star, and with a disper- 
sion corresponding to the uncertainty on the observed parallax. 
Note that this procedure is not strictly correct, as the gaussian 
distribution should in fact be centered on the (unknown) true 
parallax. To overcome this difficulty, a fully bayesian approach 
will be applied in yet another analysis of the data, as described 
next in Sect. 13. 31 

Some a priori information has nevertheless been included 
to truncate the gaussian parallax distribution in the present 
Monte Carlo approach, since the giant stars in our sample 
should not be brighter than Mnp - -2.5 (see Figs. [T3l 
I14lll6t . That threshold has been increased to Mhj, - -5 (see 
Fig|10> for the stars flagged as supergiants. This prescription 
thus corresponds to assigning a minimum admissible parallax 
TTmin to any given star It prevents that very small parallaxes 
drawn from the gaussian distribution yield unrealistically large 
space velocities. The parallax distribution used in the Monte 
Carlo simulation thus writes 



1 



(2n)i/V;,, 



■ exp 



1^ / TT - TTohs 



if n>n„ 



Pin) ^0 \f n<n„ 



(15) 
(16) 



The mean space velocities and the velocity dispersions are 
calculated for each simulated sample, and finally we adopt the 
average of these values (of the mean velocities and of the ve- 
locity dispersions) over 4000 simulated samples as the best es- 
timate of the true kinematic parameters. We thus obtain: 



{U) = -10.25 +0.15 kms-i 
(V) = -22.81 +0.15 kms"' 
{W) = -7.98 ±0.09 km s-i. 



(17) 



If there is no net radial and vertical motion in the so- 
lar neighbourhood (see however Sect. 3.3.6), we may write 
Uq - -{U) and = -(W). The results given by ( fTTt are 
very close to those estimated from stars with relative errors 
on the parallax smaller than 20% (see (|5j)- They are in agree- 
ment with the values derived by Dehnen & Binney (1998) on 
the basis of Hipparcos proper motions of main sequence stars 
(t/o = 10.00 ± 0.36kms-i, Wq = 7.17 ± 0.38kms-i). The 
value of Uq is not in perfect accordance with the one derived 
by Brosche et al. (2001), who found - 9.0+0.5 kms"' from 
photometric distances and Hipparcos proper motions of KO-5 
giants, nor with the one derived by Zhu (2000) who found t/© = 
9.6 + 0.3 kms"' with the same stars as ours but without the ra- 
dial velocity data. The value of contradicts slightly the one 



derived by Bienayme (1999), who found Wq - 6.7 + 0.2 kms"' 
from Hipparcos proper motions. Nevertheless, these consider- 
ations on the solar motion are not very useful since we stress 
in Sect. 3.3.6 that the motion of the Sun is difficult to derive 
anyway because there are conceptual uncertainties on the mean 
motion of stars in the solar neighbourhood. 
For the velocity dispersions, we find 



a-u =40.72 + 0.58 km s-i 
cry = 32.23 + 1.41 km s-i 
a-w = 22.55 + 0.95 km s-i. 



(18) 



These velocity dispersions are somewhat larger than those 
found for the restricted sample which is not surprising since 
60% of the high-velocity stars (as identified in Sect. 3.3.6) are 
not present in the sample restricted to the most precise paral- 
laxes. Regarding the asymmetric drift and the slightly larger 
age of the M giants, the Monte Carlo method and the analysis 
of the restricted sample reach the same conclusion. 

Concerning the mixed component of the velocity dispersion 
tensor in the plane, we find 



a-ly = 188.23 + 40.9 km^s-^ 



(19) 



which leads to a clearly non-zero vertex deviation of 



= 1/2 arctan(2cr^y/(cr^ - cr^)) = 16.2° + 5.6°. (20) 

Interestingly, Soubiran et al. (2003) found no vertex devia- 
tion for low-metallicity stars in the disk, and concluded that this 
was consistent with an axisymmetric Galaxy: this is absolutely 
not the case for our sample of late-type giants. Furthermore, 
the value we find for the vertex deviation is larger than the 
one derived for late-type stars by Dehnen & Binney (1998), 
who showed that the vertex deviation drops from 30° for young 
stellar populations (probably due to young groups concentrated 
near the origin of the f/V-plane, see also Sect. 3.3.7) to a con- 
stant value of 10° for older populations. Bienayme (1999) also 
found from Hipparcos proper motions a vertex deviation of 9.2° 
for the giant stars. We conclude from our sample that the vertex 
deviation is significantly non-zero: we suggest a possible ori- 
gin for this vertex deviation (which could be the same origin as 
for the vertex deviation of younger populations) in Sect. 3.3.7. 

3.3. Bayesian approach 

To obtain the kinematic characteristics of our sample in a 
more rigorous way, we have applied the Luri-Mennessier (LM) 
method described in detail by Luri (1995) and Luri et al. 
(1996). The starting point of this method is a model describing 
the basic morphological characteristics we can safely expect 
from the sample (spatial distribution, kinematics and absolute 
magnitudes), and a model of the selection criteria used to de- 
fine the sample. These models are used to build a distribution 
function (DF) intended to describe the observational character- 
istics of the sample. The a priori distribution function adopted 
is a linear combination of partial DFs, each of which describes 
a group of stars (called "base group"). Each partial DF com- 
bines a kinematic model (the velocity ellipsoid introduced by 
Schwarzschild 1907) with a gaussian magnitude distribution, 
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and an exponential height distribution uncorrelated with the 
velocities. The phenomenological model adopted is obviously 
not completely rigorous, but has the advantage of being able 
to identify and quantify the different subgroups present in the 
data and possibly related to extremely complex dynamical phe- 
nomema, which cannot be easily parametrized. 

The values of the parameters of the DF can be determined 
from the sample using a bayesian approach: the model is ad- 
justed to the sample by a maximum likelihood fit of the pa- 
rameters. The values of the parameters so obtained provide the 
best representation of the sample given the a priori models as- 
sumed. 

In the following Sects. 3.3.1 and 3.3.2 we describe the 
ingredients of the models used in this paper Thereafter, in 
Sects. 3.3.3 to 3.3.6, we present the results of the maximum 
likelihood fit, and we discuss in Sect. 3.3.7 the different possi- 
ble physical interpretations of those results. 

3.3.1. Phenomenological model 

To build a phenomenological model of our sample, we assume 
that it is a mixture of stars coming from several "base groups". 
A given group represents a fraction w, of the total sample and 
its characteristics are described by the following components: 

- Spatial distribution: An exponential disk of scale height Zq 
in the direction perpendicular to the galactic plane 



(feid, I, b) - exp 



d cosb 



(21) 



and a uniform distribution along the galactic plane (a real- 
istic approximation for samples in the solar neighbourhood 
like ours). 

Velocity distribution: A Schwarzschild ellipsoid for the ve- 
locities of the stars with respect to their Local Standard of 
Rest (LSR): 



where 

U' = {U- Uo) cos h - (V - Vo) sin Z„ 
V' ^(U - Uo) sin lr + (V- Vq) cos Z„ 



(22) 



(23) 



and where is the vertex deviation. 

Galactic rotation: An Oort-Lindblad rotation model at first 
order, where the rotation velocity is added to the mean LSR 
velocity (given by the Schwarzschild ellipsoid above). The 
same values were adopted for the Oort constants as in Sects. 
3.1 and 3.2 

Luminosity: We have adopted a gaussian distribution for the 
absolute magnitudes of the stars: 



<Pm(M) = exp 



2\ 



1/M-MoV 
2 \ o-M / ^ 



(24) 



This is just a first rough approximation. A better model 
would include a dependence of the absolute magnitude on 
a colour index but is more complex (we leave it for future 
papers). 



- Interstellar absorption: In the LM method, the coiTection 
of interstellar absorption is integrated in the formalism. A 
model, giving its value as a function of the position {d, I, b), 
is needed: the Arenou et al. (1992) model has been chosen 
here. 

The DF of each base group is simply the product of the 
space, velocity and magnitude distributions presented here, 
with different values for the model parameters. The total dis- 
tribution function D of the sample is a linear combination of 
the partial DFs of the different base groups. 

Both the relative fractions of the different base groups and 
the model parameters will be determined using a Maximum 
Likelihood fit, as described in Sect. 13.3.31 However, the num- 
ber of groups («g) composing our sample is not known a pri- 
ori. A likelihood test - like Wilk's test, Soubiran et al. (1990) 
- will be used to determine it: ML estimations are performed 



with n„ 



1,2,3, 



. and the maximum likelihoods obtained for 



each case are compared using the test to decide on the most 
Ukely value of ng. 

3.3.2. Observational selection and errors 

As pointed out above, the correct description of the observed 
characteristics of the sample requires that its selection criteria 
be included in the model. In our case, our sample of giants (like 
the Hipparcos Catalogue as a whole) is composed primarily of 
stars belonging to the Hipparcos survey plus stars added on the 
basis of several heterogeneous criteria. 

In the case of the Hipparcos survey, the selection criteria 
are just based on the apparent V magnitude. For a given line of 
sight, with galactic latitude b, the survey is complete up to 



y = 7.3-Hl.l|sinZ7|. 



(25) 



In our case, the Hp Hipparcos magnitude has been used in- 
stead of V. Therefore the survey stars in our sample will follow 
a similar "completeness law" in Hp that we have adopted to be 
approximately 



Hp = 7.5 + l.l|sin^7| 



(26) 



assuming an average Hp - V of 0.2 mag for the stars in our 
sample. 

This law, altogether with the cutoff in declination {6 > 0°), 
allows us to quite realistically model the selection of the survey 
stars in the sample, but leaves out 10% of the total sample. 
In order to be able to use the full sample, we must define the 
selection of the non-survey stars as well, as (approximately) 
done by the following completeness law: 



- Completeness linearly decreasing from 1 at Hp 
1 . 1 1 sin down to zero at Hp - 11 .5 



7.5 + 



This condition together with the completeness up to Hp = 
7.5 + 1.1| sin^l , and with the cutoff in declination, defines the 
selection of the complete sample. 

The individual errors on the astrometric and photometric 
data are also taken into account in the model. We assume that 
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the observed values are produced by gaussian distributions of 
observational errors around the "true" values, with standard de- 
viations given by the errors quoted in the Hipparcos Catalogue, 
the Tycho-2 catalogue or the CORAVEL database. 

In the end, we can define for each star a joint distribution 
function of the true and observed values: 

^i(x,z\\e)^D{x\\e)e{z\\x)S{z) (27) 

where x are the true values (position in the LSR, velocities with 
respect to the LSR, absolute magnitude), z the observed values 
(position on the sky, parallax, proper motions, radial velocity, 
apparent magnitude), the set of parameters of the model, S 
the selection function taking into account the selection criteria 
of the sample, and e the gaussian distribution of observational 
errors. 

For each star, we can then easily deduce the distribution 
function of the observed values: 

0(z||6l) = J n{x,z\\0)dx (28) 

The fact that the selection function and the error distribu- 
tion are taken into account in this distribution function prevents 
from all the possible biases. 



giants and supergiants with 'young' kinematics (further dis- 
cussed in Sect. I3.3.6t . one of high-velocity stars and finally 
one group of 'normal' stars. This third group exhibits plenty 
of small-scale structure, and this small-scale structure can be 
successfully modelled by 3 more base groups (leading to a to- 
tal of 6 groups in the sample, called groups Y, HV, HyPl, Si, 
He, and B, further discussed in Sect. 3.3.6). These supplemen- 
tary groups are statistically significant, as revealed by a Wilks 
test. Statistically, solutions with 7 or 8 groups are even bet- 
ter, but these solutions are not stable anymore and are thus not 
useful (they depend too much on the observed values of some 
individual stars). Table 1 lists the values of the ML parame- 
ters for a model using just the 5177 (out of 5397) survey stars 
which comply with the completeness law (Eq. (26)). Table 2 
lists those values for a model using the whole sample of 6030 
stars with our modified selection law. The models for the survey 
stars alone and for the whole sample give very similar results 
(except for the scale height Zq of group HV - see Sect. 3.3.6 - 
which is not well constrained since our sample does not go far 
enough above the galactic plane), thus providing a strong in- 
dication that the results for the whole sample are reliable. The 
errors on the parameters listed in Tables 1 and 2 are the disper- 
sions coming from the simulated samples. 



3.3.3. Maximum likelihood 



3.3.4. Group assignment 



The principle of maximum likelihood (ML) can be briefly de- 
scribed as follows: let z be the random variable of the observed 
values following the density law given by 0{z\\Oo), where 6q - 

(Wi , Moi , CTm, , t/Oi , cry; , Vo, , O-yj , Wo, , (TWi , Zoi , /y, , W2, . . . , lv„^ ) 

is the set of unknown parameters on which it depends. 
The Likelihood Function is defined as 



L(0) = f]0(z,||0). 



(29) 



The value of 6 which maximizes this function is the ML 
estimator, 0ml, of the parameters Bq characterizing the density 
law of the sample. It can be shown that 0ml is asymptotically 
non-biased, asymptotically gaussian and that for large samples, 
it is the most efficient estimator (see Kendall & Stuart 1979). 

Once the ML estimator of the parameters has been found, 
simulated samples (generation of random numbers following 
the given distribution) are used to check the equations and pro- 
grams developed for the estimation. They also allow a good 
estimation of the errors on the results and the detection of pos- 
sible biases: once a ML estimation 0ml has been obtained, sev- 
eral samples corresponding to the parameters 0ml are simu- 
lated and the method is applied to them. The comparison of the 
results (which could be called "the estimation of the estima- 
tion") with the original 0ml allows us to detect possible biases. 
The dispersion of these results cr(0ML) can be taken as the error 
of the estimation. 

The low maximum likelihood obtained when the full sam- 
ple is modelled with a single base group indicates that the 
kinematic properties of giant stars in the solar neighbourhood 
cannot be fitted by a single Schwarzschild ellipsoid. The first 
acceptable solution requires three base groups: one of bright 



The probabilities for a given star to belong to the various groups 
is provided by the LM method, and the most probable group 
along with the corresponding probability is listed in Table A. 1 . 
Indeed, let Wj be the a priori probability (i.e. the ML parameter 
of Table 2) that a star belongs to the yth group and 0j{z \ Oj) 
the distribution of the observed quantities z in this group (de- 
duced from the phenomenological model adopted and depend- 
ing on the parameters 0j of this model for the group). Then, 
using Bayes formula, the a posteriori probability for a star to 
belong to the jth group given its measured values z» is: 



P{* e Gj I z,) = 



Wj0(z.\0j) 

i:lUwt0t(z.\0ky 



(30) 



Using this formula the a posteriori probabilities that the star 
belongs to a given group can be compared, and the star can be 
assigned to the most hkely one (Table A.l). 

Note that this procedure, like any method of statistical clas- 
sification, will have a certain percentage of misclassifications. 
However, the reliability of each assignment is clearly indicated 
by the probability given in column 35 of Table A. 1 . 

3.3.5. Individual distance estimates 

Once a star has been assigned to a group, and given the 
ML estimator of the group parameters and the observed values 
for the star z*, one can obtain the marginal probability density 
law "Rid) for the distance of the star from the global probability 
density function. 
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Table 1. Maximum-Likelihood parameters obtained using 5177 survey stars only. The number of stars belonging to each group 
(as listed in column 34 of Table A.l) is given in the last row. Although obtained from the assignment process described in 
Sect. 13.3.11 it is not, strictly, a ML-parameter, in contrast to the 'fraction' w, {i = 1, ...6) of the whole sample belonging to the 
considered group, as listed in the previous row. Therefore, the observed fraction of stars in each group could be slightly different 
from the true fraction w, of that group in the entire population. 





group Y 


group HV 


group HyPl 


group Si 


group He 


group B 


Mo (mag) 


0.61 + 0.91 


2.21 ± 0.22 


1.15 ± 0.37 


1.08 ± 0.57 


1.27 ±0.37 




1.23 ± 0.17 


(Tm (mag) 


I.OJ ± 0.24 


1 in 1 n 1 n 
i./9 ± U.IO 


n OO 1 ATI 

0.06 ± O.zl 


1.16 ± 0.24 


0.88 ±0.12 




1 nc 1 n no 
l.UD ± O.Oo 


Uo (km s-') 


-11.63 ±2.30 


-17.56 ± 3.46 


-31.23 ±0.99 


5.19 ±3.03 


-42.09 ± 4.76 




-2.92 ± 1.54 


(Tt// (km s"') 


15.61 ± 1.22 


53.45 ±2.14 


11.46 ± 1.53 


13.48 ± 2.66 


26.08 ± 13.13 




31.82 ± 1.56 


Vo(kms-i) 


-11.57 ± 2.08 


-43.96 ± 3.20 


-20.00 ± 0.75 


4.22 ± 1.95 


-50.60 ± 4.67 




-15.15 ±2.36 


<Tv' (km s"') 


9.44 + 1.36 


36.14 ± 1.79 


4.89 ± 1.10 


4.59 ± 4.22 


8.62 ± 3.03 




17.60 ± 0.78 


Wo (km s-') 


-7.79 ± 0.95 


-7.76 ± 2.30 


-4.82 ± 1.85 


-5.57 ±2.13 


-6.94 ± 3.56 




-8.16 ±0.91 


(Tn/ (km s"') 


6.94 ± 0.78 


32.54 ± 1.54 


8.76 ± 1.08 


9.43 ± 2.23 


16.66 ± 1.62 




16.26 ± 0.77 


Zo (pc) 


73.44 ± 6.33 


901.04 ±440.13 


128.70 ± 19.21 


149.43 ± 23.67 


201.80 ±55.24 


196.05 ± 12.34 


/v n 


17.07 ± 5.81 


0.22 ± 4.89 


-5.65 ± 4.59 


-14.17 ±25.28 


-5.67 ± 11.16 




-0.22 ± 5.21 


fraction w 


0.086 ± 0.013 


0.149 ±0.022 


0.071 ± 0.007 


0.050 ± 0.030 


0.065 ± 0.017 




0.579 ± 0.047 


number of stars 


345 


505 


334 


204 


372 




3417 


Table 2. Same as Table 1 for the full sample of 6030 stars 




group Y 


group HV 


group HyPl 


group Si 


group He 


group B 


Mo (mag) 


0.68 ± 0.25 


2.04 ± 0.23 


0.95 ± 0.13 


0.87 ± 0.19 


1.21 ± 0.18 


1.04 ± 0.07 


a-M (mag) 


1.78 ±0.06 


1.43 ± 0.08 


0.74 ± 0.10 


1.12 ± 0.08 


0.95 ± 0.09 


1.07 ± 0.03 


Uo (km s-') 


-10.41 ±0.94 


-18.50 ±2.82 


-30.34 ± 1.54 


6.53 ± 1.93 


-42.13 ± 1.95 




2.78 ± 1.07 


(Ty> (km s"') 


15.37 ± 1.09 


58.02 ± 1.87 


11.83 ± 1.34 


14.37 ± 2.05 


28.35 ± 1.68 


33.30 ± 0.70 


Vo(km s-i) 


-12.37 ±0.90 


-53.30 ±3.10 


-20.27 ± 0.62 


3.96 ± 0.67 


-51.64 ± 1.07 




15.42 ± 0.82 


CTv' (km s"') 


9.93 ± 0.75 


41.36 ± 1.70 


5.08 ± 0.76 


4.63 ± 0.71 


9.31 ± 1.22 


17.94 ± 0.80 


Wo (km s-') 


-7.75 ± 0.57 


-6.61 ± 1.82 


-4.82 ± 0.80 


-5.80 ± 1.15 


-8.06 ± 1.30 




8.26 ± 0.38 


cTw (km s"') 


6.69 ± 0.58 


39.14 ± 1.69 


8.75 ± 0.74 


9.66 ± 0.82 


17.10 ± 1.63 


17.65 ± 0.34 


Zo (pc) 


80.27 ±6.18 


208.08 ± 27.89 


106.27 ± 14.37 


127.96 ± 19.65 


132.89 ± 9.42 


141.21 ± 5.20 


iv n 


16.40 ± 10.26 


0.16 ±4.77 


-8.79 ± 4.05 


-11.96 ± 3.53 


-6.53 ± 2.77 




2.18 ± 1.60 


fraction w 


0.096 ± 0.008 


0.106 ±0.010 


0.070 ± 0.008 


0.053 ± 0.009 


0.079 ± 0.009 


0.596 ± 0.015 


number of stars 


413 


401 


392 


268 


529 


4027 



This can then be used to obtain the expected value of the dis- 
tance 

d'R(d)dd (31) 



and its dispersion 

6^ = I {d- df nd) dd. (32) 
Jo 

The first can be used as a distance estimator free from biases 
and the second as its eiTor (Columns 28 and 29 of Table A.l). 
Fig - El reveals that the biases in the distance derived from the 
inverse parallax are a combination of those resulting from trun- 
cations in parallax and apparent magnitude as discussed in Luri 
& Arenou(1997). 

3.3.6. The various groups present in tine sample 

We have plotted in Fig.|5]the stars in the f/V-plane by using 
the expected values of U and V deduced from the LM method. 
The 6 different groups are represented on this figure. The struc- 
ture of the f/V-plane is similar to the one already identified (in 
Fig.0 for the stars with precise parallaxes io-„/n < 20%). This 



similarity is a strong indication that the subgroups identified by 

the LM method are not artifacts of this method. Several groups 

can be identified with known kinematic features of the solar 

neighbourhood. 

Group Y: The young giants 

The first group is one with 'young' kinematics (yellow group 
on Fig. |9}, with a small velocity dispersion and with a vertex 
deviation of U, - 16.4°. Quite remarkably, the value of (Mnp) 
for this group is 0.68 (see Table 2), the brightest among all 6 
groups. An important property of this group is thus its young 
kinematics coupled with its large average luminosity: the most 
luminous giants and supergiants ' in our sample thus appear 
to have a small dispersion in the t/V-plane, in agreement with 
the general idea that younger, more massive giants are predom- 
inantly found at larger luminosities in the Hertzsprung-Russell 
diagram and are, at the same time, concentrated near the origin 
of the t/V-plane. It must be stressed that nothing in the LM 
method can induce such a correlation between velocities and 
luminosity artificially. This result must thus be considered as 
a robust result, even more so since the group Y of young gi- 
ants was present in all solutions, irrespective of the number of 

' almost all stars a priori flagged as supergiants (see Sects. 2.2, 3.2 
and Fig. 5) indeed belong to group Y 
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Fig. 8. Comparison of the (biased) distances obtained from a 
simple inversion of the parallax, and the maximum-likelihood 
distances d obtained from the LM method. 



groups imposed. This group is centered on the usual antisolar 
motion (Dehnen & Binney 1998) in U and W as seen in Table 2. 

Fig. Allocates stars from group Y in the HR diagram, con- 
structed from the V - I colour index as provided by the colour 
transformation of Platais et al. (2003), based on the measured 
Hp - Vt2 colour index (and is thus more accurate than the 
Hipparcos V - I index; see Appendix for more details). The 
isochrones of Lejeune & Schaerer (2001), for a typical metal- 
licity Z = 0.008, indicate that the age of stars from group Y is 
on the order of several 10^ to a few 10** yr. 

It must be remarked at this point that the observed mem- 
bers of group Y, as displayed on Fig. ^| appear to be 
brighter than the true average absolute magnitude for the 
group listed in Tables 1 and 2 ({Mhj,} = 0.68) and as- 
signed by the LM method. This is a natural consequence of 
the Malmquist bias (Malmquist 1936) for a magnitude-limited 
sample. Nevertheless, even the true average absolute magni- 
tude makes group Y the brightest and youngest one. 
Group HV: The high-velocity stars 

This group is composed of high-velocity stars (represented 
in blue on Fig. Those stars are probably mostly halo or 
thick-disk stars, even though the value of the scale height of this 
population is poorly constrained (hence the large difference be- 
tween that parameter derived from the full sample - Table 2- 
and from survey stars only - Table 1) because our sampling 
distance is too small: indeed Fig. |S] shows that the distance of 
most stars in our sample is smaller than the scale height of the 
thick disk (665 pc < Zo tMck < 1000 pc). This group represents 
about 10% of the whole sample (Table 2, although, like Zq, this 
parameter is not well constrained), and this value is consistent 
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Fig. 9. All the stars plotted in the f/V-plane with their values of 
U and V deduced from the LM method. The 6 different groups 
are represented on this figure: group Y in yellow, group HV in 
blue, group HyPl in red, group Si in magenta, group He in green 
and group B in black. Note that the yellow group (Y) extends 
just far enough to touch both the red (HyPl) and magenta (Si) 
groups. 



with the mass fraction of the thick disk relative to the thin disk 
in the solar neighbourhood. 

Going back to the different results obtained for the veloc- 
ity dispersions in Eqs. (8) and (18), notice that 247 stars (out 
of 401) of this group have relative errors on their Hipparcos 
parallax higher than 20%. 

Clearly we see on Fig. ^2 that stars of this group are old, 
mostly older than 1 Gyr and it seems clear that some stars at 
the bottom of the HR diagram are as old as the Galaxy itself. 
It is also striking that the envelope of this kinematically 'hot' 
group seems to correspond to a portion of circle approximately 
centered on a zero-velocity frame with respect to the galactic 
center, and with a radius of the order of 280 km s ' . We refrain, 
however, from providing numerical values for the position of 
this circular envelope in the f/V-plane, because of the rather 
limited number of stars defining this limit and the absence of 
antirotating halo stars in our sample. Note that a lower limit to 
the local escape velocity is provided by the velocity (f/^ H- (V -H 
220)2 + Wyi^ = 281 km s"' of HIP 89298, the fastest star in 
group HV. 

The detection of this group of high-velocity stars 'cleans' 
the sample and allows us to study the fine structure of the ve- 
locity distribution of disk stars. 
Group HyPl: Hyades-Pleiades supercluster 

The group represented in red in Fig. |9] occupies the well- 
known region of the Hyades and Pleiades superclusters. The 
large spatial dispersion of the HyPl stars clearly hints at their 
supercluster rather than cluster nature, since they are spread aU 
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Fig. 10. HR diagram of group Y. The absolute magnitude used 
here is My. Isochrones of Lejeune & Schaerer (2001) for Z = 
0.008 and log(age(yr))= 7.4, 7.6, 8, 8.3, 8.55, 8.75, 8.85, 9. 
V -I indices were computed from the colour transformation of 
Platais et al. (2003) 




V-l 



Fig. 11. HR diagram of group HV. Isochrones of Lejeune & 
Schaerer (2001) for Z = 0.008 and log(age(yr))= 8.3, 8.55, 
8.75, 8.85, 9, 9.3, 9.45. 9.5, 9.6, 9.7. V -I indices were com- 
puted from the colour transformation of Platais et al. (2003) 

over the sky with a wide range of distances (up to 500 pc). 
The Hyades ({U} = -40 km s"', (V) = -25 km s"', see 
Dehnen 1998) and Pleiades ({U) = -15 km s"', (V) = 
-25 km s"', see Dehnen 1998) superclusters are known since 
Eggen (1958, 1975). It was also noticed (Eggen 1983) that sev- 
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Fig. 12. HR diagram of group HyPl. Isochrones of Lejeune 
& Schaerer (2001) for Z = 0.02 and log(age(yr))= 8.3, 8.55, 
8.75, 8.85, 9, 9.3, 9.45. 9.5, 9.6, 9.7. We calculated V-I using 
the colour transformation of Platais et al. (2003). The stars at 
the bottom-right of the diagram could be very old metal rich 
stars (in contradiction with the age-metallicity relation which is 
obviously not correct if no other factors are taken into account). 



eral young clusters (NGC 2516, IC 2602, a Persei) have the 
same V-component as those superclusters. We obtain (V) - 
-20.3 km s-ifor our supercluster structure, and thus refine the 
old value (V) = -25 km s"\ with in addition the fact that the 
group is tilted in the anti-diagonal direction of the t/V-plane, 
i.e., it has a negative vertex deviation (Z,, - -8.7°). The concept 
of branches in the f/V-plane (Skuljan et al. 1999), of quasi- 
constant V but slightly tilted in the anti-diagonal direction, will 
be further discussed below in relation with the other groups. 

Metallicity is available for 17 stars of the HyPl group 
(Mc William 1990), and it is interesting to remark that the 
average metallicity seems to be close to solar ([Fe/H]= 0; 
Fig.ll5> as already noticed by Chereul & Grenon (2001). This 
is quite unusual for giant stars. The metallicity of the stars from 
group B for example is centered on [Fe/H]= -0.2. Plotting 
the HyPl stars in a HR diagram (Fig. ll2> and using the solar- 
metallicity isochrones of Lejeune & Schaerer (2001) reveals 
that the stars forming the HyPl group are by far not coeval (de- 
spite the fact that HyPl stars have a smaller age spread than the 
Si and He groups discussed below, especially because the HyPl 
group is lacking young supergiants). Although a precise deter- 
mination of the ages of HyPl stars would require the knowledge 
of their metallicities, it seems nevertheless clear that there are 
a lot of clump stars with ages reaching 1 Gyr, and that there are 
clearly stars older than 2 Gyr, in sharp contrast with the ages 
of about 80 and 600 Myr for the Pleiades and Hyades clusters 
themselves. This result of age heterogeneity, already noticed 
by Chereul 8c Grenon (2001) for the Hyades supercluster, is an 
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Fig. 13. HR diagram of group Si. Isochrones of Lejeune & 
Schaerer (2001) for Z = 0.008 and log(age(yr))= 8.3, 8.55, 
8.75, 8.85, 9, 9.3, 9.45. 9.5, 9.6, 9.7. 



important clue to identify the origin of supercluster-like struc- 
tures, as discussed in Sect. 3.3.7. 

Finally, the value (W) - -4.8 + 0.8 km s"' differs signifi- 
cantly from the vertical solar motion {Wq = 7 to 8 km s"'; see 
Sects. 13.21 and group B in Tables 1 and 2), indicating that 
the group has a slight net vertical motion. 
Group Si: Sirius moving group 

The Sirius moving group (represented in magenta in Fig.|5) 
is known since Eggen (1958, 1960) and is traditionally located 
at (U) = 10 km s"', (V) = -5 km s"' (Dehnen 1998). The LM 
method refines those values and locates it at (U) = 6.5 km s ', 
(V) - 3.9 km s The spatial distribution once again indi- 
cates that this group has a supercluster-like nature. Metallicity 
from Mc William (1990) is available for 12 stars of the Si group 
(Fig. I15> . Contrarily to the situation prevailing for the HyPl 
group, the metallicity distribution within the Si group appears 
similar to that for the bulk of the giants in the solar neighbour- 
hood, as represented by group B. Isochrones for a typical value 
of Z = 0.008 (Fig.ll3> indicate that the ages are widely spread, 
even more so than for the HyPl group. 
Group He: Hercules stream 

The 'Hercules stream' (represented in green in Fig.|5), located 
by the LM method at {U) - -42 km s"^ (V) = -51 km s^'has 
been named by Raboud et al. (1998) the M-anomaly. It corre- 
sponds to a global outward radial motion of the stars which 
lag behind the galactic rotation: known since Blaauw (1970), 
it is traditionally centered on {U) - -35 km s"', (V) = 
-45 km s"'(see Fux 2001). It is strongly believed since several 
years that its origin is of dynamical nature, as further discussed 
in the next section. Once again, the Hertzsprung-Russell dia- 
gram of Fig. [3] indicates a wide range of ages in this group 
(for a typical value of Z = 0.008). 
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Fig. 14. HR diagram of group He. Isochrones of Lejeune & 
Schaerer (2001) for Z = 0.008 and log(age(yr))= 8.3, 8.55, 
8.75, 8.85, 9, 9.3, 9.45. 9.5, 9.6, 9.7. 



Note that the three groups HyPl, Si and He could be in- 
terpreted in term of extended branches crossing the f/V-plane 
(Skuljan et al. 1999, Nordstrom et al. 2004): one is Hercules, 
another is the combination of Hyades and Pleiades, and the 
third one is Sirius. The high dispersion of the streams azimuthal 
(V) component confirms this view. The branches are somewhat 
tilted along the anti-diagonal direction. These properties could 
be understood in the classical theory of moving groups as dis- 
cussed in next section but also in the context of a dynamical 
origin for the substructure. 
Group B: Smooth Background 

Most stars are part of an 'axisymmetric', 'smooth' background 
represented in black in Fig. |9l The average metallicity of this 
group seems to be slightly subsolar (Fig. \\5\ as expected for 
a sample of disk giants: Girardi & Salaris (2001) found an 
average metallicity of [Fe/H]=-0.12 + 0.18, while for their 
sample of F and G dwarfs, Nordstrom et al. (2004) found 
[Fe/H]=-0.14 + 0.19. Using isochrones for a typical value of 
Z = 0.008 (see Fig.ll6>. we see that there is a large spread in 
age. This is typical of the mixed population of the galactic disk, 
composed of stars born at many different epochs since the birth 
of the Galaxy. 

In the t/V-plane, the velocity ellipsoid of this group is not 
centered on the value commonly accepted for the antisolar mo- 
tion: it is centered instead on {U) - -2.78 ± 1.07 km s~^ 
However, the full data set (including the various superclusters) 
does yield the usual value for the solar motion (see Sect. 3.2). 
This discrepancy clearly raises the essential question of how to 
derive the solar motion in the presence of streams in the so- 
lar neighbourhood: does there exist in the solar neighbourhood 
a subset of stars having no net radial motion? If the smooth 
background is indeed an axisymmetric background with no net 
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Fig. 15. Histogram of the metallicity in groups B, HyPl and Si 
for the stars present in the analysis of Mc William ( 1990) 




Fig. 16. HR diagram of group B. Isochrones of Lejeune & 
Schaerer (2001) for Z = 0.008 and log(age(yr))= 8.3, 8.55, 
8.75, 8.85, 9, 9.3, 9.45. 9.5, 9.6, 9.7. 

radial motion, we have found a totally different value for the so- 
lar motion. Nevertheless, we have no strong argument to assess 
that this is the case, especially if the 'superclusters' have a dy- 
namical origin, as proposed in the next section. Moreover, the 
group Y of young giants is centered on the commonly accepted 
value of (U) - -10.41 + 0.94 km s"' and this difference be- 
tween groups Y and B clearly prevents us from deriving with- 
out ambiguity the solar motion. If the giant molecular clouds 
(GMC's) from which the young stars arose are on a circular or- 



bit, then the (U) value of group Y is the acceptable one for the 
antisolar motion, but nothing proves that the GMC's are not 
locally moving outward in the Galaxy under the effect of the 
spiral pattern. No value can thus at the present time be given 
for the radial solar motion but only some different estimates 
depending on the theoretical hypothesis we make on the nature 
of the substructures observed in velocity space. Theoretical in- 
vestigations and dynamical simulations thus appear to be the 
only ways to solve this problem. On the other hand, the value 
of (W) = -8. 26+0. 38km s"' forgroupBis in accordance with 
the usual motion of the sun perpendicular to the galactic disk 
(see Sect. 3.2), and thus seems to be a reliable value because 
streams have a smaller effect on the vertical motion of stars. 



3.3.7. Physical interpretation of tine groups: dynamical 
streams? 

Several mechanisms may be responsible for the substruc- 
trure observed in velocity space in the solar neighbourhood. 
Hereafter, we list them and confront them with the results of 
our kinematic study. 

A first class of mechanisms is that associated with inhomo- 
geneous star formation responsible for a deviation from equi- 
librium in the solar neighbourhood: this theory states that a 
large number of stars are formed (almost) simultaneously in a 
certain region of the Galaxy and create a cluster-like structure 
with a well-defined position and velocity. After several galac- 
tic rotations, the cluster will evaporate and form a tube called 
'supercluster'. Stars in the 'supercluster' still share common V 
velocities when located in the same region of the tube (for ex- 
ample in the solar neighbourhood) for the following reason, 
first put forward by Woolley (1961): if the present galactocen- 
tric radius of a star on a quasi-circular epicycUc orbit equals 
that of the sun (denoted -Wq), and if such a star is observed with 
a peculiar velocity v - V + Vq, then its guiding-center radius 
TUg writes 



2B 



(33) 



where Xg is the position of its guiding-center in the cartesian 
reference frame associated with the LSR (in the solar neigh- 
bourhood approximation, the impact of yg is negligible), and B 
is the second Oort constant. Woolley (1961) pointed out that 
disk stars (most of which move on quasi-circular epicyclic or- 
bits) which formed at the same place and time, and which 
stayed together in the Galaxy after a few galactic rotations 
(since they are all currently observed in the solar neighbour- 
hood) must necessarily have the same period of revolution 
around the Galactic center, and thus the same guiding-center 
m^, and thus the same velocity V - v - Vq according to Eq. 
i33\ . This theory has thus the great advantage of predicting 
extended horizontal branches crossing the f/V-plane (similar 
to those observed in our sample on Fig. |5] and already identi- 
fied by Skuljan et al. 1999, and by Nordstrom et al. 2004 in 
their sample of F and G dwarfs). Moreover, it is easy to under- 
stand in this framework that the Group HyPl seems to be more 
metal-rich than the smooth background (see Fig. ll5l l. However, 
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this theory does not explain the tilt of the branches that we ob- 
serve in our sample. Moreover, to explain the wide range of 
ages observed in those branches (see Figs. ^] and also 
Chereul & Grenon 2001), the stars must have formed at dif- 
ferent epochs (see e.g. Weidemann et al. 1992) out of only two 
large molecular clouds (one associated with the Sirius branch 
and another with the Hyades-Pleiades branch). Chereul et al. 
(1998) suggested that the supercluster-like velocity structure 
is just a chance juxtaposition of several cluster remnants, but 
this hypothesis requires extraordinary long survival times for 
the oldest clusters (with ages > 2 Gyr) in the supercluster-like 
structure. This long survival time is made unlikely by the argu- 
ment of Boutloukos & Lamers (2003) who found that clusters 
within 1 kpc from the Sun having a mass of m x IO'^Mq can 
survive up to m Gyr in the Galaxy: indeed heavy clusters of 
more than 2x10"* Mq are probably very rare in the disk. 

An explanation for the wide range of ages observed in the 
superclusters could be the capture of some older stars by the 
high concentration of mass in a molecular cloud, at the time of 
the formation of a new group of stars (while other stars would 
be scattered by the molecular cloud) . Though this hypothesis 
could be tested in A^-body simulations involving gas, it would 
imply that the local clumps of the potential not only perturb 
the motion of stars but dominate it, which is in contradiction 
with the high predominance of the global potential on the local 
clumps. This explanation is thus highly unlikely too. 

The second class of mechanisms involves the theory of hi- 
erarchic formation of the galaxies; following this theory, galax- 
ies were built up by the merging of smaller precursor struc- 
tures. It is known since the discovery of the absorption of the 
Sagittarius dwarf galaxy by the Milky Way (Ibata et al. 1994) 
that some streams in the Galaxy are remnants of a merger with 
a satellite galaxy. Helmi et al. (1999) showed that some debris 
streams are also present in the galactic halo near the position 
of the Sun. Such streams could also be present in the velocity 
substructure of the disk: it seems to be the case of the Arcturus 
group at t/ ^ km s"', V - -115 km s"', as recently pro- 
posed by Navarro et al. (2004). In this scenario, the streams 
observed in our kinematic study could also be the remnants 
of merger events between our Galaxy and a satellite galaxy. 
The merger would have triggered star formation whereas the 
oldest giant stars would be stars accreted from the companion 
galaxy. A merger with a satellite galaxy would moreover in- 
duce a perturbation in W, as we observe in group HyPl in which 
{W} = -4.8 +0.8 km s"'. If the hierarchical ('bottom-up') cos- 
mological model is correct, the Milky Way system should have 
accreted and subsequently tidally destroyed approximately 100 
low-mass galaxies in the past 12 Gyr (see Bullock & Johnston 
2004), which leads to one merger every 120 Myr, but the 
chance that two of them (leading to the Hyades-Pleiades and 
Sirius superclusters) have left such important signatures in the 
disk near the position of the Sun is statistically unlikely (al- 
though not impossible). 

The third class of mechanisms is the class of purely dy- 
namical ones. A dynamical mechanism that could cause sub- 
structure in the local velocity distribution is the disturbing ef- 
fect of a non-axisymmetric component of the gravitational po- 
tential, like the rotating galactic bar. The Hercules stream was 



recently identified with the bimodal character of the local ve- 
locity distribution (Dehnen 1999, 2000) due to the rotation of 
the bar if the Sun is located at the bar's outer Lindblad reso- 
nance (OLR). Indeed, stars in the Galaxy will have their orbits 
elongated along or perpendicular to the major axis of the bar 
(orbits respectively called the LSR and OLR modes in the ter- 
minology of Dehnen 2000), depending upon their position rel- 
ative to the resonances, and both types of orbits coexist at the 
OLR radius. Moreover, all orbits are regular in a 2-dimensional 
(2D) axisymmetric potential, but the perturbation of the triax- 
ial bar will induce some chaos. Fux (2001) showed that in the 
region of the Hercules stream in velocity space, the chaotic re- 
gions, decoupled from the regular regions, are more heavily 
crowded. The Hercules stream has thus also been interpreted 
as an overdensity of chaotic orbits (Fux 2001) due to the ro- 
tating bar Quillen (2003) has confirmed that when the effect 
of the spiral structure is added to that of the bar, the 'chaotic' 
Hercules stream remains a strong feature of the local DF, whose 
boundaries are refined by the spiral structure. This stream in 
the local DF seems thus related to a non-axisymmetric pertur- 
bation rather than to a deviation from equilibrium due to in- 
homogeneous star formation. Following Miihlbauer & Dehnen 
(2003), the Galactic bar naturally induces a non-zero vertex de- 
viation on the order of 10° and the vertex deviation found in 
Sects. 3.1 and 3.2 would thus be partly related to the Hercules 
stream. Nevertheless, the other streams present in the data are 
also partly responsible for the vertex deviation. The likely dy- 
namical origin of the Hercules stream has been the first ex- 
ample of a non-axisymmetric origin for a stream in velocity 
space: the other streams could thus be related to other non- 
axisymmetric eff'ects. Chereul & Grenon (2001) proposed that 
the Hyades supercluster represents in fact an extension of the 
Hercules stream, but the LM method has shown that the two 
features are clearly separated in the f/V-plane. 

De Simone et al. (2004) have shown that the structure of the 
local DF could be due to a lumpy potential related to the pres- 
ence of transient spiral waves in the shearing sheet (i.e. a small 
portion of an infinitesimally thin disk that can be associated 
with the solar neighbourhood). Julian & Toomre (1966) have 
studied, in the shearing sheet, the response of a stellar disk to a 
point-like perturbation of the potential. They found that this re- 
sponse is in the form of density waves whose wavecrests swing 
around from leading to trailing. These waves are transient be- 
cause their amplitudes are amplified transitorily and then die 
away. The duration of these transient spiral waves is of the or- 
der of an orbital period, i.e. of the order of 10*^ years, and they 
impart momentum to stars: this mechanism tends to put stars in 
some specific regions of the f/ V-plane in the simulations of De 
Simone et al. (2004), thus creating streams as observed in our 
sample. These spiral waves cause radial migration in the galac- 
tic disk near their corotation radius, while not increasing the 
random motions and preserving the overall angular momentum 
distribution (Sellwood & Binney 2002). The seemingly pecu- 
liar chemical composition of the group HyPl (i.e., a metallic- 
ity higher than average for field giants, as suggested by the 
17 giant stars analyzed by Mc William 1990 and displayed in 
Fig. El ^Iso reported by Chereul & Grenon 2001) thus sug- 
gests that the group has a common galactocentric origin in the 
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inner Galaxy (where the interstellar medium is more metal-rich 
than in the solar neighbourhood) and that it was perturbed by a 
spiral wave at a certain moment. This specific scenario (Pont et 
al., in preparation) would explain why this group is composed 
of stars sharing a common metallicity but not a common age. 
The group Si could also be a clump recently formed by the 
passage of a transient spiral. 

The simulations of De Simone et al. (2004) can create 
streams with a range of 3 Gyr or more in age, and this is thus 
a mechanism that can explain that main result of our study. 
Another characteristic of the simulations is that they tend to re- 
produce the observed branches, and their origin ultimately lies 
in the same mechanism as that elucidated by Woolley (1961). 
However, the tilt of the branches in the t/V-plane is not repro- 
duced by the simulations. 

A question which arises in the framework of this dynami- 
cal scenario is the following: is it by chance that such a large 
number of young clusters and associations are situated on the 
same branches in the f/V-plane (more than half of the OB as- 
sociations listed in Table A.l of de Zeeuw 1999 are situated in 
the regions of the HyPl and Si groups on Fig.O? Probably not, 
which suggests that the same transient spiral which gave their 
peculiar velocity to the HyPl and Si groups has put these young 
clusters and associations along the same UV branches. For ex- 
ample, it is quite striking that the Hyades cluster itself is metal- 
rich with a mean [Fe/H]= 0.13 (Boesgaard 1989), thus pointing 
towards the same galactocentric origin as the HyPl group. The 
entire cluster could have been shifted in radius while remain- 
ing bound since the effect of a spiral wave on stars depends on 
the stars' phase with respect to the spiral, and the phase does 
not vary much across the cluster. Moreover, since most clus- 
ters and associations are young, they should not have crossed 
many transients, which suggests that one and the same spiral 
transient could have formed some clusters and associations (by 
boosting star formation in the gas cloud) and could at the same 
time have given them their peculiar velocity in the f/V-plane. 
The relation between spiral waves and star formation is indeed 
well established (e.g. Hernandez et al. 2000 who found a star 
formation rate S FR{t) with an oscillatory component of period 
500 Myr related to the spiral pattern). As a corollary, we con- 
clude that the dynamical streams observed among K and M 
giants are young kinematic features: integrating backwards (in 
a smooth stationnary axisymmetric potential) the orbits of the 
stars belonging to the streams makes thus absolutely no sense, 
and reconstructing the history of the local disk from the present 
data of stars in the solar neighbourhood becomes tricky. In this 
dynamical scenario, the deviation from dynamical equilibrium 
that is present among samples of young stars is closely related 
to the deviation from axisymmetry existing in the Galaxy. Of 
course, we do not exclude that the position of some clusters 
and OB associations in the same region of velocity space as the 
dynamical streams could be the result of chance (Chereul et al. 
1998). 

If this dynamical scenario is coiTect, the term 'dynamical 
stream' for the branches in velocity space seems more appro- 
priate than the term 'supercluster' since they are not caused by 
contemporaneous star formation but rather involve stars that 



do not share a common place of birth: stars in the streams just 
share at present time a common velocity vector 

It should be noted that those non-axisymmetric perturba- 
tions, as well as the minor mergers, could lead to some asym- 
metries in the spatial distribution of stars in the galactic disk 
on a large scale (see Parker et al. 2004). Since our sample does 
not cover the whole sky, and is anyway restricted to the solar 
neighbourhood, we are not in a position to detect those asym- 
metries. 

To conclude this section, let us stress that the dynamical, 
non-axisymmetric theory to explain the substructure observed 
in the velocity space is largely preferred over the theory which 
views it as remnants of clusters of stars sharing a common ini- 
tial origin, essentially because of the wide range of ages of the 
stars composing the identified subgroups. In fact, we stress that 
the deviation from equilibrium among young stars (Dehnen & 
Binney 1998) and the dynamical streams exhibited in this pa- 
per could be closely related to each other. The presence of dy- 
namical streams in our sample of K and M giants is clearly 
responsible for the vertex deviation found for late-type giants, 
but we even suggest that the same origin could hold as well 
for early-type stars. It is indeed quite striking on Fig.|9lthat the 
group Y extends just far enough in the UV -plane to touch both 
the HyPl and Si branches. The presence of stars from these 
two streams in group Y (sent in that region of the f/V-plane 
at the time of their formation because of the peculiar velocity 
imparted by the spiral wave that created them) imposes a very 
specific value to the vertex deviation (ranging from 15° to 30° 
and more), in agreement with the high values often observed 
for young stars (see Dehnen & Binney 1998). This idea that 
the vertex deviation for younger populations could in fact have 
the same dynamical origin as the vertex deviation for old ones 
was already proposed by Mayor (1972, 1974). Here, we also 
argue that even the specific initial conditions of young groups 
of stars could be due to the same phenomenon. 

However, nature is of course not so simple and the fea- 
tures of the DF are presumably related to a mixture of several 
phenomena. Notably, the initial conditions in the simulations 
should be more complex than the simple 2D Schwarzschild ve- 
locity ellipsoid used by De Simone et al. (2004): in fact, purely 
axisymmetric substructure could already be present in the solar 
neighbourhood (see the structure of the f/V-plane in Dejonghe 
& Van Caelenberg 1999). Other phenomena that could have 
an influence on the structure of velocity space are the follow- 
ing: a triaxial or clumpy dark halo, giant molecular clouds, and 
close encounters with the Magellanic Clouds (Rocha-Pinto et 
al. 2000). Theoretical investigations in this area should thus 
clearly be pursued. 

4. Summary and perspectives 

This paper presented the first kinematic analysis of 53 1 1 K and 
719 M giants (after excluding the binaries for which no center- 
of-mass velocity could be estimated) in the solar neighbour- 
hood which included radial velocity data from an important 
survey performed with the CORAVEL spectrovelocimeter We 
also used proper motions from the Tycho-2 catalogue. 
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First, we analyzed the kinematics of the sample restricted 
to the 2774 stars with parallaxes accurate to better than 20%, 
and then we made full use of the 6030 available stars and evalu- 
ated the kinematic parameters with a Monte-Carlo method. We 
found that the asymmetric drift is larger for M giants than for 
K giants due to the fact that the M giants must be a little older 
than the K giants on average. We also found the usual value for 
the solar motion when assuming that the whole sample has no 
net radial and vertical motion (Sect. 3.2). 

Then a maximum-likelihood method, based on a bayesian 
approach (Luri et al. 1996, LM method), has been applied to the 
data and allowed us to derive simultaneously maximum likeli- 
hood estimators of luminosity and kinematic parameters, and 
to identify subgroups present in the sample. Several subgroups 
can be identified with known kinematic features of the solar 
neighbourhood (namely the Hyades-Pleiades supercluster, the 
Sirius moving group and the Hercules stream). Isochrones in 
the Hertzsprung-Russell diagram reveal a very wide range of 
ages for stars belonging to these subgroups. We thus conclude 
that this substructure in velocity space is probably related to the 
dynamical perturbation by transient spiral waves (as recently 
modelled by De Simone et al. 2004) rather than to cluster rem- 
nants. The emerging picture is thus one of dynamical streams 
pervading the solar neighbourhood. A possible explanation for 
the presence of young group/clusters in the same area of the 
f/V-plane as those streams is that they have been put there by 
the spiral wave associated with their formation, while kinemat- 
ics of the older stars of our sample have also been disturbed 
by the same wave. The seemingly peculiar chemical composi- 
tion of the Hyades-Pleiades stream (also reported by Chereul & 
Grenon 2001) suggests that this stream originates from a spe- 
cific galactocentric distance and that it was perturbed by a spi- 
ral wave at a certain moment and radially pushed by the wave in 
the solar neighbourhood (see Sellwood & Binney 2002; Pont et 
al. in preparation). This would explain why this stream is com- 
posed of stars sharing a common metallicity but not a common 
age. A careful metallicity analysis of this stream would be of 
great interest in order to confirm this scenario. The Sirius mov- 
ing group could also be a feature recently formed by the pas- 
sage of a transient spiral, while the Hercules stream would be 
related to the bar's outer Lindblad resonance (Dehnen 1999, 
2000, Fux 2001). The position of all these streams in the UV- 
plane is responsible for the vertex deviation of 16.2° + 5.6° for 
the whole sample (see Sect. 3.2). We even argue that the ver- 
tex deviation observed among large samples of early-type stars 
(see Dehnen & Binney 1998) and the specific kinematic initial 
conditions of some young open clusters and OB associations 
could in fact have the same dynamical origin as those streams 
of giants. A better understanding of the streams should start 
with a chemical analysis of the stars belonging to them. As a 
first step, their photometric indices could be investigated. Note 
also that an important consequence of the dynamical origin of 
the streams is that it makes no sense to integrate backwards (in 
a smooth stationary axisymmetric potential) the orbits of the 
stars belonging to them. 

The group of background stars (group B) has a distribution 
in velocity space close to a Schwarzschild ellipsoid but is not 
centered on the classical value found for -t/o in Sect. 3.2 when 



considering the full sample, including the streams Instead we 
find(f/) = -2.78+ 1.07km s"'. This discrepancy clearly raises 
the essential question of how to derive the solar motion in the 
presence of dynamical perturbations altering the kinematics of 
the solar neighbourhood (the net radial motion of stars in the 
solar neighbourhood can be of the order of 10 km s ' in the 
simulation of De Simone et al. 2004): does there exist in the 
solar neighbourhood a subset of stars having no net radial mo- 
tion which can be used as a reference against which to measure 
the solar motion? 

Theoretical investigations in this area should thus clearly be 
pursued, and in particular dynamical modeling. We have shown 
that the fine structure of phase space in the solar neighbour- 
hood cannot be interpreted in terms of an axisymmetric steady- 
state model. Nevertheless, an axisymmetric model revealing 
all the fine structure of the axisymmetric distribution func- 
tion in the solar neighbourhood (Dejonghe & Van Caelenberg 
1999, Famaey et al. 2002, Famaey & Dejonghe 2003) would 
give ideal initial conditions (more complex than a simple 2D 
Schwarzschild velocity ellipsoid; see especially the structure 
of the f/V-plane in Dejonghe & Van Caelenberg 1999) for 3D 
A^-body simulations that could afterwards reproduce some non- 
axisymmetric features observed in the solar neighbourhood 
(see Fux 1997, De Simone et al. 2004 for 2D simulations). 
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values are found in columns 9 to 12; the kinematic study 
made then use of the Hipparcos proper motions instead of 
the Tycho-2 ones. In most cases, the absence of a Tycho-2 
proper motion is caused by the fact that the star has a close 
visual companion or is too bright for the Tycho detection. 
13-16 Same as columns 9 to 12 for Hipparcos proper motions. 
17 The normalized absolute difference between the Hipparcos 
and Tycho-2 proper motions: A/i = lyUeip - A'Tyc2l/e//, where 



+ e; 



Vtk-: 



_^)'^^, where 



yu = (jjIcos'^6 + yU^)'/2 and 
e, denotes the standard error of quantity /. The quantity Afi 
may be used as a diagnostic tool to identify long-period bi- 
naries (Kaplan & Makarov 2003). The basic idea behind 
this tool is the following. For binary stars with orbital pe- 
riods much longer than the duration of the Hipparcos mis- 
sion, the proper motion recorded by Hipparcos is in fact 
the vector addition of the true proper motion and of the or- 
bital motion. This orbital motion averages out in the Tycho- 
2 value, since it is derived from measurements spanning a 



Table A. 1 is only available in electronic form at the CDS 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 



http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ /. 
first page of Table A. 1 is printed hereafter. 



Nevertheless, the 
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much longer time base (on the order of a century, as com- 
pared to 3 yr for Hipparcos). Therefore, a difference be- 
tween the Tycho-2 and Hipparcos proper motions (beyond 
the combined error bar encapsulated by e^) very likely hints 
at the binary nature of the star This diagnostic will be fully 
exploited in a forthcoming paper devoted to the binary stars 
present in our sample. Note that Fig. 1 presents the his- 
togram of AyU, which indeed reveals the presence of an ab- 
normally large tail at A// > 1 .5. 25- 
18-19 Hipparcos Hp magnitude and associated standard error 

20 Tycho-2 Vti magnitude. 

21 Hp - Vt2- For visual binaries (flag 4 in column 24), this 
colour index has not been listed and the value 0.0 is given 
instead, because the Hp magnitude appears to be a compos- 
ite value for the two visual components. In most of these 
cases, the Vt2 magnitude of the visual companion is given 
in the last column of the table. For large amplitude variables- 
stars, the Hp - Vji index is meaningless as well, and as been 
set to 0. 

22 y - / colour index in Cousins' system as provided by 
field H40 of the Hipparcos catalogue. This index is use- 
ful for constructing the Hertzsprung-Russell (HR) diagranBl- 
of the sample. It should be stressed that the H40 field of the 
Hipparcos catalogue does not provide a directly measured 
quantity. It has instead been computed from various colour 
transformations based on the B - V index from field H3734- 
which neither is a directly measured quantity. 

23 y - / colour index in Cousins' system as provided by the 
colour transformation based on the measured Hp - Vti 
colour index (Platais et al. 2003). This value is thus in prin- 
ciple more reliable than the Hipparcos H40 value (note, 
however, that the colour transformation provided by Platais 
et al. (2003) had to be extrapolated somewhat, since it is 
provided in the range -2.5 < Hp - Vti ^ -0.20, whereas 
our data set goes up to Hp - Vti = 0.1). It has been used 
to draw the HR diagram of the sample. For large-amplitude 
variable stars, the median V - I index has been taken di- 
rectly from Platais et al. (2003), instead of being computed 
from the colour transformation based on Hp - Vti (set to 0. 
in those cases). 

24 Binarity flag: 

*: no evidence for radial- velocity variations; 
0: spectroscopic binary (SB), with no orbit available. The 
star had to be discarded from the kinematic analysis; 
1 : SB with an orbit available, or with a center-of-mass ve- 
locity which can be estimated reasonably well. Column 25 
then contains the system's center-of-mass velocity. The ref- 
erence with the orbit used to derive the center-of-mass ve- 
locity is listed in the last column. If no reference is given, 
the center-of-mass velocity has been estimated from the 
available CORAVEL data; 

2: supergiant star, with a substantial radial-velocity jitter 
(see Fig. 5); 

3: uncertain case: either SB or supergiant; 

4: visual binary with a companion less than 6" away (as 

listed by the Tycho-2 catalogue); 

5,6,7,8: as 0,1,2,3 but for a visual binary. It should be 
stressed here that visual binaries have not been searched for 



exhaustively among our target stars. Whenever the Tycho- 
2 catalogue lists a companion star less than 6" away from 
the target star, the binarity flag has been set to 4. In most 
of these cases, the Hipparcos Hp magnitude corresponds to 
the composite magnitude and the Tycho-2 proper motions 
are identical for the two components. The Vt:2 magnitude 
of the companion is listed in the last column of the table; 
9: binary supergiant with an orbit available. 

26 Average radial velocity (based on CORAVEL observations) 
or center-of-mass velocity for SBs (the last column pro- 
vides the bibliographic code of the reference providing the 
orbit used), and standard error (set to 0.3 km s"' in the case 
of center-of-mass velocity). 

27 The absolute Hp magnitude, coiTected for interstellar red- 
dening (according to the model of Arenou et al. 1992), and 
based on the LM distance listed in column 28. 

29 The maximum-likelihood distance based on the LM esti- 
mator (see Sect. 13. 3> . and its associated standard error (see 
Eq. (|33). 

30 The interstellar absorption Ay, based on the LM distance 
and the model of Arenou et al. (1992). 

33 The U, V and W components of the heliocentric space 
velocity deduced from the LM method (corrected for 
the galactic diff'erential rotation to first order using A - 
14.82 kms-'kpc"' and B = -12.37kms-ikpc"'). 

35 The most likely group to which the star belongs, and the as- 
sociated probability. The various groups are the following: 

1 or Y: stars with a young kinematics; 

2 or HV: high- velocity stars; 

3 or B: stars defining the smooth background (the most pop- 
ulated group); 

4 or HyPl: stars belonging to the Hyades-Pleiades stream; 

5 or Si: stars belonging to the Sirius stream; 

6 or He : stars belonging to the Hercules stream. 

36 In the case of a spectroscopic binary with an avail- 
able orbit: Bibliographic code (according to the standard 
ADS/CDS coding) of the reference listing the orbit used 
as the source of the center-of-mass velocity (The code 
2005A&A...???..???J refers to the forthcoming paper by 
Jorissen et al. 2005 devoted to the analysis of the binary 
content of the present sample). In the case of a visual bi- 
nary, the Vj2 magnitude of the companion. An asterisk in 
that column indicates the presence of a note in the remark 
file. 

It must be stressed that columns 27 to 35 contain model- 
dependent data, as they were derived by the LM method. They 
depend upon, e.g., the particular choice for the values of the 
Oort constants, the interstellar extinction model, the a priori 
choice of the various distribution functions,... 
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Table A.l. The left-hand side of the first page of Table A.l, fully available in electronic form at the CDS. See text for a description 
of the column contents. 
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-23.4 
-11.6 
-18.7 
-44.5 
18.1 
5.6 
-7.9 
57.8 
-0.3 
99.1 
9.9 
-8.1 
18.7 
0.0 
95.1 
1094 
16.0 
7.8 
34.7 
-7.0 
82.7 
83.3 
11.9 
131.5 
35.9 
17.7 
8.4 
103.5 
130.9 
0.0 
-6.3 
4.3 
-29.2 
-6.6 
-24.7 
0.1 
0.0 
-37.3 
-10.8 
73.3 
59.5 
87.1 
17.5 
-13.7 
17.9 
19.8 



1.2 
1.0 
1.1 
1.3 
0.9 
1.1 
0.8 
0.8 
1.2 
1.2 
0.9 
1.0 
1.1 
1.0 
1.3 
1.1 
0.0 
1.2 
1.1 
0.8 
1.1 
0.8 
1.1 
1.1 
1.0 
0.9 
1.2 
1.4 
0.9 
1.0 
1.1 
1.3 
0.0 
1.3 
0.9 
1.3 
1.2 
1.2 
1.0 
0.0 
1.0 
0.8 
0.9 
1.5 
0.9 
1.5 
0.9 
1.1 
0.7 



-0.4 
1.3 
15.6 
18.0 
-17.0 
-33.6 
3.6 
-29.0 
13.1 
-12.0 
-52.8 
-29.4 
-46.7 
-6.4 
-11.8 
0.9 
0.0 
24.7 
-13.1 
-0.4 
-1.8 
-16.8 
-4.1 
-18.3 
8.2 
2.1 
-119.4 
-7.8 
-33.5 
-6.4 
55.1 
37.0 
0.0 
-6.0 
10.6 
-4.1 
-16.0 
-0.6 
1.9 
0.0 
11.3 
-8.3 
-13.3 
10.8 
3.8 
324 
-25.3 
-4.5 
9.4 



1.3 
1.1 
1.1 
1.3 
1.0 
1.1 
0.9 
0.8 
1.3 
1.3 
0.9 
1.0 
1.1 
1.1 
1.3 
1.1 
0.0 
1.3 
1.1 
0.8 
1.2 
0.8 
1.1 
1.0 
1.1 
0.9 
1.2 
14 
0.9 
1.1 
1.0 
1.4 
0.0 
1.2 
0.9 
1.4 
1.1 
1.2 
1.1 
0.0 
1.1 
0.9 
1.1 
1.5 
0.9 
1.5 
1.1 
1.2 
0.7 



61.89 
-4.88 
49.71 
-20.78 
-10.61 
-18.88 
-44.73 
18.50 
4.74 
-8.91 
58.45 
0.04 
98.95 
7.86 
-8.01 
20.68 
-15.80 
95.58 
112.32 
16.36 
8.82 
34.93 
-5.09 
85.60 
82.91 
12.29 
131.20 
35.24 
18.87 
7.09 
104.61 
128.28 
-2.50 
-6.09 
4.16 
-30.17 
-5.04 
-25.63 
0.04 
-1.90 
-37.27 
-10.76 
73.95 
60.54 
86.82 
18.32 
-13.46 
16.16 
19.02 



0.84 
1.49 
0.85 
0.67 
0.60 
0.66 
0.57 
0.71 
0.65 
0.50 
0.74 
0.77 
0.66 
0.56 
0.63 
1.07 
1.51 
0.48 
0.71 
0.67 
0.63 
0.61 
0.63 
0.68 
0.74 
0.95 
0.83 
0.71 
0.75 
0.73 
0.83 
0.69 
1.80 
0.98 
0.81 
0.84 
0.81 
0.75 
0.94 
1.87 
0.53 
0.85 
0.65 
0.46 
0.66 
0.78 
0.65 
0.60 
0.82 



-0.22 
-0.20 
14.27 
18.19 

-18.37 

-34.10 
3.74 

-28.01 
12.31 

-11.05 

-52.87 

-31.13 

-46.75 
-8.15 

-10.94 
1.42 

-30.38 
23.80 

-14.52 
-0.82 
-1.82 

-17.05 
-4.69 

-16.84 
10.65 
1.16 
-117.65 
-8.78 

-33.40 
-7.99 
54.77 
38.30 
2.02 
-5.73 
9.10 
-4.08 

-16.73 
-0.52 
3.36 
-0.81 
11.00 
-6.87 

-15.28 
10.39 
3.73 
33.73 

-24.35 
-2.94 
8.66 



0.56 
0.69 
0.60 
0.55 
0.51 
0.46 
0.53 
0.47 
0.53 
0.52 
0.53 
0.59 
0.61 
0.46 
0.60 
0.81 
1.09 
0.45 
0.47 
0.57 
0.68 
0.49 
0.49 
0.52 
0.58 
0.74 
0.52 
0.55 
0.60 
0.56 
0.51 
0.62 
2.03 
0.68 
0.58 
0.76 
0.54 
0.51 
0.67 
1.82 
0.55 
0.65 
0.53 
0.45 
0.46 
0.69 
0.56 
0.56 
0.53 



-0.30 
1.09 
-0.10 
-0.94 
041 
0.30 
0.17 
-0.44 
-0.54 
-0.09 
0.34 
1.01 
-0.06 
-0.28 
-0.37 
0.98 
0.00 
0.13 
1.73 
0.26 
0.53 
0.23 
-0.68 
1.48 
-0.06 
0.15 
-0.71 
-0.19 
0.30 
0.07 
0.47 
-1.01 
0.00 
-0.16 
-0.89 
0.43 
0.09 
0.48 
0.77 
0.00 
-0.07 
-0.53 
0.62 
0.43 
-0.19 
0.66 
-0.57 
-1.11 
-0.73 
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Table A. 1. The right-hand side of the first page of Table A. 1 , fully available in electronic form at the CDS . For the sake of clarity, 
the HEP number has been repeated in the leftmost colunm of the printed version, and the notes have been Usted at the bottom of 
the table. 
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Hp 
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Vt2 
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new 


B 
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(km/s) 


Mhp 


d ± 
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V 

(km/s) 


W 
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D 


P(g) 
(%) 


21 


7.6923 


0.0012 


7.683 


0.009 


1.43 


1.24 


* 


-11.72 


0.16 


1.29 


184.4 


28.0 


0.07 


-42.5 


-31.4 


-1.6 


3 


49.7 


31 


7.7648 


0.0020 


7.793 


-0.028 


1.46 


1.47 


* 


3.52 


0.21 


-0.67 


471.6 


131.2 


0.07 


6.5 


6.2 


-0.6 


5 


61.6 


36 


7.8503 


0.0015 


7.805 


0.045 


1.11 


1.00 





0.00 


0.00 


0.00 


0.0 


0.0 


0.00 


0.0 


0.0 
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0. 


73 


7.0364 


0.00 11 


7.054 


-0.018 


1.44 


1.41 


* 


-11.91 


0.20 


-0.51 


245.6 


35.3 


0.60 


24.9 


-1.6 


24.9 


3 


90.4 


119 


6.9834 


0.0010 


6.934 


0.049 


1.04 


0.97 


* 


0.53 


0.21 


0.95 


146.4 


18.1 


0.20 


10.3 


1.9 


-9.6 


5 


58.7 


121 


6.9957 


0.0011 


6.964 


0.032 


1.16 


1.09 


* 


-15.96 


0.21 


-0.33 


272.1 


48.0 


0.15 


43.0 


-20.2 


-21.9 


3 


97.7 


123 


7.3790 


0.0013 


7.342 


0.037 


1.13 


1.05 


* 


-13.32 


0.22 


0.36 


203.0 


26.3 


0.49 


42.0 


7.3 


9.3 


3 


90.6 


136 


6.7635 


0.0010 


6.722 


0.041 


1.12 


1.02 




-7.69 


0.21 


0.31 


182.6 


25.0 


0.15 


-5.6 


-24.8 


-17.7 


3 


77.9 


142 


7.4726 


0.0013 


7.404 


0.069 


0.78 


0.83 




-21.68 


0.24 


3.29 


63.7 


2.9 


0.16 


6.9 


-20.4 


2.1 


3 


68.3 


181 


6.7960 


0.0034 


6.980 


-0.184 


2.46 


2.30 


* 


-21.39 


0.23 


-1.98 


377.9 


67.0 


0.89 


25.8 


-11.4 


-17.8 


1 


54.4 


258 


6.7616 


0.0008 


6.713 


0.049 


1.08 


0.97 


* 


-18.50 


0.21 


0.56 


162.5 


16.3 


0.15 


-17.5 


-56.3 


-23.8 


6 


71.4 


282 


8.0651 


0.0018 


8.075 


-0.010 


1.39 


1.36 


* 


-47.28 


0.22 


-0.73 


535.3 


144.5 


0.15 


47.3 


-73.3 


-21.8 


2 


44.8 


288 


8 0688 


0.0013 


8.028 


0.041 


1.04 


1.03 


* 


3.33 


0.38 
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32.8 


0.23 
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-49.2 


-58.6 


2 


49.9 
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6.3544 


0.0045 
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2.28 


2.22 
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0.00 


0.00 


0.0 


0.0 


0.00 


0.0 


0.0 
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0. 


323 


7.5413 


0.0010 


7.475 


0.066 


0.98 


0.85 


* 


-16.05 


0.39 


0.92 


192.9 


34.6 


0.19 


11.5 


-11.8 


-5.6 


3 


76.8 


368 


8.0579 


0.0016 


8.045 


0.013 


1.21 


1.21 


* 


-16.06 


0.22 


0.83 


270.3 


67.3 


0.07 


-19.9 


-18.6 


9.9 


3 


52.0 


374 


7.1831 


0.0011 


7.255 


-0.072 


0.96 


1.73 


4 


-23.98 


0.29 


0.27 


216.8 


56.2 


0.23 


29.3 


-22.9 


-13.0 


3 


97.4 


379 


5.8511 


0.0006 


5.794 


0.057 


1.02 


0.92 


* 


-28.43 


0.20 


0.61 


99.5 


5.7 


0.25 


-27.1 


-47.2 


1.0 


6 


64.3 


399 


6.3941 


0.0010 


6.415 


-0.021 


1.34 


1.43 





0.00 


0.00 


0.00 


0.0 


0.0 


0.00 


0.0 


0.0 


0.0 





0. 


466 


7.3999 


0.0009 


7.358 


0.042 


1.14 


1.02 


* 


-9.72 


0.22 


-0.39 


314.7 


51.6 


0.30 


-14.6 


-20.3 


-3.3 


4 


46.2 


472 


7.1505 


0.0009 


7.115 


0.035 


1.15 


1.06 





0.00 


0.00 


0.00 


0.0 


0.0 


0.00 


0.0 


0.0 


0.0 





0. 


496 


7.0181 


0.0012 


7.067 


-0.049 


1.58 


1.60 


* 


-14.36 


0.22 


0.37 


199.1 


27.1 


0.15 


-18.3 


-32.9 


-9.4 


3 


57.3 


502 


6.7361 


0.0007 


6.774 


-0.038 


1.58 


1.53 


* 


9.98 


0.20 


-1.14 


350.5 


77.3 


0.15 


5.7 


9.9 


-9.0 


5 


59.9 


504 


7.6749 


0.0012 


7.608 


0.067 


1.07 


0.85 


* 


-10.40 


0.21 


1.15 


188.2 


20.7 


0.15 


-54.5 


-49.0 


-19.4 


6 


82.7 


525 


7.0996 


0008 


7.062 


0.038 


1.05 


1.05 


* 


-73.82 


0.25 


1.49 


121.7 


11.1 


0.18 


-14.1 


-83.7 


22.6 


2 


94.9 


598 


8.1267 


0.0015 


8.126 


0.001 


1.44 


1.29 


* 


29.49 


0.23 


0.09 


378.5 


109.4 


0.15 


-20.4 


12.5 


-20.3 


3 


92.4 


607 


7.4919 


0.0011 


7.438 


0.054 


0.99 


0.94 


* 


16.44 


0.22 


2.14 


110.4 


10.0 


0.14 


-41.5 


-47.7 


-71.0 


2 


90.8 


625 


7.5541 


0.0008 


7.503 


0.051 


0.98 


0.96 


* 


-15.02 


0.51 


0.96 


188.0 


20.8 


0.23 


-16.3 


-29.0 


-6.1 


4 


49.3 


704 


8.0043 


0.0013 


7.947 


0.057 


1.09 


0.91 


* 


-37.60 


0.22 


1.05 


237.8 


38.2 


0.07 


11.2 


-56.0 


-0.4 


6 


46.4 


714 


6.9335 


0.0012 


6.990 


-0.056 


1.59 


1.64 


* 


-11.78 


0.54 


-1.14 


370.0 


82.0 


0.24 


3.8 


-19.4 


-8.3 


1 


57.5 


716 


6.4099 


0.0006 


6.354 


0.056 


0.97 


0.92 


* 


16.11 


0.20 


0.71 


129.3 


12.6 


0.15 


-72.5 


-1.4 


6.6 


3 


92.4 


755 


7.1836 


0.0010 


7.128 


0.056 


0.93 


0.93 


* 


5.00 


0.22 


3.15 


59.7 


2.4 


0.15 


-37.0 


-14.0 


4.5 


3 


63.9 


779 


9.4783 


0.0041 


9.750 


-0.272 


2.22 


2.68 


9 


-74.31 


0.13 


-5.78 


3601.3 


771.2 


2.47 


158.0 


-51.3 


40.7 


1 


59.2 


799 


8.5304 


0.0016 


8.471 


0.059 


1.04 


0.90 


* 


-1.34 


0.53 


1.03 


306.6 


100.3 


0.07 


11.2 


-2.1 


-3.1 


3 


51.6 


824 


7.8382 


0.0022 


7.985 


-0.147 


2.15 


2.13 


* 


18.91 


0.22 


-0.08 


357.1 


87.1 


0.15 


-13.7 


19.4 


1.6 


3 


86.5 


834 


8.2793 


0.0043 


8.263 


0.016 


1.00 


1.19 


1 


-38.93 


0.30 


0.67 


251.4 


47.7 


0.61 


46.7 


-17.2 


-0.7 


3 


98.2 


852 


7.0964 


0.0048 


7.345 


-0.249 


2.64 


2.59 


* 


-14.17 


0.34 


-0.53 


312.1 


72.7 


0.15 


19.4 


-16.4 


-10.1 


3 


86.2 


858 


6.9629 


0.0007 


6.970 


-0.007 


1.42 


1.34 


* 


-15.05 


0.26 


-0.05 


226.9 


37.5 


0.24 


28.8 


0.5 


10.9 


3 


64.2 


868 


8.3296 


0.0013 


8.298 


0.032 


1.13 


1.09 


* 


-40 89 
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498 'i 


196 7 


0.07 


9.0 


-18.4 


36.9 


3 


80.2 


871 


9.1126 


0.0393 


9.133 


-0.020 


0.73 


1.42 


6 


-6.50 


0.30 


-1.39 


934.2 


319.8 


0.65 


17.1 


-1.4 


-1.3 


1 


38.0 




6.8469 


0.0010 






1.17 


1.13 




-16.90 


0.31 


-0.42 


221.4 


30.7 


n ^4 


40 fi 


^ 4 


1 Q 9 


3 


89.8 


909 


7.8542 


0.0016 


7.806 


0.048 


1.12 


0.98 


* 


-24.72 


0.22 


0.29 


315.9 


89.6 


0.07 


17.8 


-14.4 


15.3 


3 


97.3 


918 


7.3180 


0.0012 


7.264 


0.054 


0.98 


0.94 


* 


0.32 


0.21 


1.05 


169.9 


18.0 


0.12 


-53.6 


-24.9 


-19.4 


3 


61.1 


967 


6.3113 


0.0009 


6.336 


-0.025 


1.48 


1.45 


* 


14.25 


0.20 


-0.07 


174.1 


20.1 


0.17 


-51.0 


-9.4 


-2.4 


3 


74.4 


969 


6.7432 


0.0006 


6.688 


0.055 


0.97 


0.93 




8.57 


0.21 


0.73 


148.5 


15.6 


0.15 


-54.1 


-21.2 


-12.6 


3 


56.8 


972 


8.3410 


0.0017 


8.297 


0.044 


1.02 


1.01 
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0.00 


0.00 


0.0 


0.0 


0.00 


0.0 


0.0 


0.0 





0. 


980 


7.5299 


0.0011 


7.484 


0.046 


1.04 


0.99 


* 


4.36 


0.14 


1.31 


162.6 


21.6 


0.16 


12.7 


5.4 


-18.1 


5 


51.5 


982 


7.7598 


0.0015 


7.758 


0.002 


1.28 


1.28 


* 


-19.64 


0.21 


-0.17 


333.4 


75.6 


0.31 


-14.4 


-31.5 


-6.3 


3 


48.1 


989 


7.7254 


0.0065 


7.919 


-0.194 


2.39 


2.35 


* 


-10.26 


0.59 


-0.17 


354.7 


77.7 


0.15 


-33.9 


-15.0 


13.4 


3 


71.4 



Notes: 

779 Very distant star (first-order correction to the differential galactic rotation inaccurate); binary of VV Cep type (Mleplb + B2.6V) 

834 2002A&A...395..885D 

871 12.154,1988A&A...207...37A 



